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- r ■ Abstract 

(N. 

t-H . We introduce a useful and rather simple class of BKP tau functions which which we shall call "easy 

tau functions" . We consider two versions of BKP hierarchy, one we will call "small BKP hierarchy" 

04 \ (sBKP) related to O(oo) introduced in [4, and "large BKP hierarchy" (1BKP) related to 0(2oo + 1) 

(— ( . introduced in |13| (which is closely related to the large O(2oo) DKP hierarchy (1DKP) introduced 

in [B]). Actually "easy tau functions" of the sBKP hierarchy were already considered in [55], here we 
are more interested in the 1BKP case and also the mixed small-large BKP tau functions [13] . Tau 
, ■ functions under consideration are equal to certain sums over partitions and to certain multi-integrals 

£N) 1 over cone domains. In this way they may be applicable in models of random partitions and models 

of random matrices. Here is the first part of the paper where sums of Schur and projective Schur 
functions over partitions are considered. 

Q-i 
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S: 

1 Introduction 

■ In the seminal papers of Kyoto school KP hierarchies related to different symmetry groups were intro- 

duced. In such a way DKP and BKP hierarchies appeared as KP hierarchies related to Z?oo and 
type root systems, while the original KP hierarchy of Dryuma-Zakharov-Shabat was assigned to the root 
system of type A^. However different realizations of these hierarchies are possible. The BKP and DKP 
hierarchies presented in [5] are subhierarchies of the standard KP one: it is related to 0(oo) subgroup 
of GL(oo) symmetry group of KP. Authors of 13] refer these BKP and DKP hierarchies as respectively 
neutral BKP and DKP hierarchies. We shall call them respectively small BKP (sBKP) and small DKP 
(sDKP) hierarchies. There is also different DKP hierarchy presented in the paper [jy Q and related to 
• • ■ (9(2oo) D GL(oo) which contains KP as the subhierarchy. We shall this hierarchy as large DKP one 

\ (1DKP). At last in [13 the BKP hierarchy related to O(2oo + 1) D GL(oo) which also contains KP as 

a subhierarchy was introduced. We shall refer it as large BKP hierachyj These " large" hierarchies are 
rather interesting and much less studied than the "small" ones. In [5] and |13) the fermionic represen- 
tation for sDKP, sBKP, 1DKP and 1BKP tau function was written down and bilinear equation ('Hirota 
equations') were presented. The tau function of these hierarchies appeared in a number of various prob- 
lems. In the paper [T3] it was shown that under certain restrictions 1BKP (and 1DKP) tau functions 
coincide with the so-called Pfaff lattice tau function [3J which in particular describes the orthogonal 
ensemble of random matrices. In |14) nice fermionic representations for orthogonal and for symplectic 
ensembles were found and in this way it was shown that the partition function of these ensembles are 
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examples of 1BKP tau function. In the recent paper [53] the coupled "large" 2-DKP hierarchy was 
introduced and the quasi-classical limit of the 1DKP hierarchy and of the 2-1DKP hierarchy was studied. 

General solutions of 1DKP and 1BKP hierarchy may be found as solutions to Hirota-type equations 
[6] , |13pl : Hirota equations for 2-1BKP hierarchy were written down in [24] . 

In the present paper we shall study certain 'simple' classes of solutions of the 1BKP and 2-1BKP 
hierarchies ( "easy tau functions" ) singled out by equations (I199P , (I200p , (I20ip as it is explained in Section[S] 
Actually these tau functions on the one hand generalize two examples presented by J. van de Leur in [T4"] 
on the other hand generalize tau functions considered in [29] and called tau functions of hypergeometric 
type. We believe that such tau functions will have various applications. In addition we find its natural 
to consider certain solutions of the 1BKP hierarchy coupled to sBKP one. Let us mark that special 
solutions of sBKP were studied in [TUJ, [TT], [3T]. sBKP was used in studies of various random processes, 
see [BT], [B7J, [B3], [BJ. 

A 1BKP tau function depends on the same set of higher times t = (ti, . . . ) as a KP tau function 
and on two discrete parameters (discrete times) I and V (instead of one parameter in the KP case) , and 
may be written in form of Schur function expansion 

Tw {t) = a *(t)n A M (1) 

AGP 

where I') are certain Pfaffians. In (JTJ sum runs over the set of all partitions denoted by P. In case 
of 2-1BKP hierarchy a typical series has the following form 

r«,(t,t)= s A (t)n AM (/,/')s M (t) (2) 

which is an analog of the Takasaki series for 2D TL hierarchy, which is 

7l(t,t)= *A(t)7M0*|.(t) (3) 

where are certain determinants. 

In the present paper we shall derive (JIJ from the fcrmionic representation of the 1BKP tau function 
given in [13) and consider a set of examples and applications. In particular we will introduce a certain 
class of 1BKP tau functions which may be considered as a generalization of the hypergeometric function 
(compare to [3D] and [SI] ) which depends on the 1BKP higher times t = (tx, <2, . . . ) and a set of parameters 
denoted by U — {U m , m G Z}, 

r(t) = £ e-^*A(t) 

ASP 

and more general tau functions, see Section[5l Here s\ are the Schur functions, t are 1BKP higher times. 
The sum ranges over the set of all partitions denoted by p. 
An example of such tau functions is 

r(t)= s ^ 

AC (m N ) 

where the sum ranges over all Young diagrams A which can be arranged into a m by TV rectangular where 
m, N are given numbers. At first sight one can think that it is just a particular case of the well-known 
series for the KP tau functions [7J, [BJ 

A 

This guess is not right: for simplicity take m = N = oo; then in KP case the numbers w\ should solve 
Plucker relations while tt\ all equal to 1 do not solve. 

3 There is an interesting remark by J.Harnad that the Hirota equation for 1DKP and 1BKP may be treated as analogues 
of the Plucker relations for isotropic Grassmannians called Cartan relations |15| . see |16| for more details on the topic of 
Cartan relations. 
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Another interesting example of the 1BKP tau function is as follows. Consider a subset of all 
partitions P denoted FP ("fat partitions") which consists of all partitions of even length of form 
(Ai, Ai, A 2 , A 2 , A 3 , A 3 , . . . ) = A U A. Then 

r(t) = ^saW 
agfp 

which is also an example of the 1BKP tau function and which we will relate to a discrete version of j3 = 4 
ensembles. This tau function will be used in Section [7] in a problem related to random motion. 

We also present different examples of tau functions which are written as multiple integrals over a cone 
domain. Such integrals appear in the theory of random matrices [17] , |89] . Let us point out the pioneer 
paper [5] , [3] which relates the orthogonal ensemble to Pfaff lattice and also the very helpful paper by J/ 
van de Leur who have shown that both ensembles of real (3—1 and f3 — 4 are examples of 1BKP theory. 
In [T] we modify some results of [H] to the case of sBKP and considered also the cases of three different 
p = 1, f3 = 2 , and /3 = 4 ensembles 

We shall explain what is the meaning of "independent variables" U — {Ui} and what are equations 
with respect to these variables. It is suitable to parameterize U by new variables t*, see (|105[) and (|106[) , 
then for certain specifications of t (see ([M]) - (|102p ) we find that 

r(t,U(t*)) = J2 e ' Uxit ' )s ^) 
Aep 

are again 1BKP tau functions now with respect to parameters t* which play the role of higher times. Let 
us mark that this tau function turned out to be a partition function for a model of random turn motion 
of vicious walkers introduced by M.Fisher [70] , This problem is considered in Section [JJ 

These times may be also considered as group times for convolution flows [22] and related to the action 
of vertex operators. Hamiltonians of these flows act in a diagonal way on the basis of Schur functions 0- 

These convolution flows on arbitrary 1BKP tau function may be also interpreted in terms of 'dual' 
multisoliton 1BKP tau functions whose higher times t* are related to parameters U = U(t*) mentioned 
above (see Section [S]). This link between two 1BKP tau functions is quite similar to the case studied in 
KP where such link between two tau functions was used for technical purposes in papers [UJ and [S5] 
and was clarified in [27] and in [25] . 

We found it is pertinent to present certain small BKP tau functions such as 

r(t',U)= £ e' Ux Q x {\t') 
agdp 

and also 1BKP tau functions coupled to sBKP tau functions (Section ref), the simplest example 

r(t,t')= £ SA(t)Qx-(|t') 

l(X)<N 

where Q\- are projective Schur functions, A~ denotes a strict partition whose parts are shifted parts of 
a partition A. This expression is a 1BKP tau function with respect to the set t = (tx,t%, • • • ) of higher 
times. At the same time it is sBKP tau function with respect to the times t'. 

2 Sums of Schur functions 

Subsets of partitions. In the following, we consider sums over partitions and strict partitions , which 
will be denoted by Greek letters a, (3. Recall [IB] that a strict partition a is a set of integers (parts) 
(a%, . . . , Qife) with ai > ■ • • > «fc > 0. The length of a partition a, denoted 1(a), is the number of 
non- vanishing parts, thus it is either k or k — 1. 

Let P be the set of all partitions. We shall need two special subsets of P. 

The first one consists of all partitions A = (Ai, . . . , A 2n ), < n e Z, A 2 ,i > 0, which satisfy 

Ai + A 2n+ i_i is independent of i , i = 1, . . . , 2n , 

4 First similar Hamiltonians were considered in the study of generalized Kontsevich model in 1401 . 
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or equivalently 

hi + hm+\-i = 2c is independent of i (hence = hi + ft-2n > 2n — 1), £ = 1, . . . , 2n , (4) 

where /ii = Ai — i + 2n , and 2c is a natural number conditioned by 2c > 2rt. This subset consists of 
all partitions A of length /(A) < 2n whose Young diagram satisfies the property that its complement 
in the rectangular Young diagram Y corresponding to (Ai + \2n) 2n coincides with itself rotated 180 
degrees around the center of Y. This set of partitions will be denoted by SCP(c) or simply SCP, for 
"self-complementary partitions" . If we introduce 

hi + ft-2n-l /rN 

yi~hi-c, c= - , (5) 

then relation ^ may be rewritten as 

Vi + V2n+l-i = . (6) 

The second subset we need consists of the partitions A which satisfy, equivalently, 

^21 = A 2l -i , i = 1,2, . . ., (7) 

or A = fi U a :— (/ii, /ii, fi 2 , M2, •• • , MfeiMfe) = M2, • • • , A*fe) 6P), or that the conjugate partitions 
of A are even, i.e., the ones whose parts are even numbers. This set of partitions will be denoted by FP, 
for "fat partitions". 

Following |TB] we will denote by DP the set of all strict partitions (partitions with distinct parts), 
namely, partitions (a±, a 2 , ■ ■ ■ , ctk), 1 < k G Z with the strict inequalities a.\ > a 2 > • • • > a*, > 0. 
Strict partitions a with the property 

a<2i = a2i-i + 1 for 2i - 1 < 1(a) , (8) 

where we set ct2i = if 1(a) = 2i — 1, will be called fat strict partitions. The set of all fat strict partitions 
will be denoted by FDPO 

The set of all self-complementary strict partitions will be denoted by SCDP. 

Let Rnm denote the set of all partitions whose Young diagram may be placed into the rectangle 
TV x M, namely, Rnm is the set of all partitions A restricted by the conditions Ai < M and £(X) < N. 



Sums over partitions. Consider the following sums (for t := (^1,^3,...), t* := (t*,^,...), t := 
(it, h, ■ ■ ■ ), t := (t u t 3 , . . . ), t* := (tl, t§, . . . ), t := (t x , i 3 , ...)), N). 

SW(t,N;U,A) := ^ A hW e~ u w s x (t) (9) 

AGP 

e(x)<N 

where h(X) = Aj — i + N. The factors Ah on the right-hand side of ([9]) are determined in terms a pair 
(A, a) =: A where A is an infinite skew symmetric matrix and a an infinite vector. For a strict partition 
h = (hi, ■ ■ ■ , ft at), the numbers Ah are defined as the Pfaffian of an antisymmetric 2n x 2n matrix A as 
follows: 

A h := Pf[l] (10) 

where for N = 2n even 

Ay = -Iji := A huhj , 1 < i < j < 2n (11) 

and for N = 2n — 1 odd 

,1, .1,, : if " -" ' (12) 

I a>hi if 1 < i < i = 2n. 

In addition we set Aq = 1. 
Then 



3 This subset was used in |26l where it was denoted by DP'. 
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where U„ , n — 0, 1, 2, . . . is a set of given complex numbers. This set is denoted by U. 

As we see the factor e" u ^ can be included into the factor Ah by redefinition of the data A as follows: 



A _v A p -u n -u m , -U n 



However we prefer to keep U as a set of parameters. 

Example We choose the following matrix A is given by 

(sgn(i-fc) if 1 < i, fc < i fl iffc<£ 

A lk = (A ) ik := { a k = \ . (14) 

[0 otherwise [0 otherwise 

Remark 1. The matrix Ai is infinite. However if in series (J9]l we put U„ = +oo for n > L, it will be the same 
as if we deals with the finite L by L matrix A, given by (|14p . 

Example 1 

Aik = (Ai) ik := 1, i<k, a k = 1 (15) 

Then 

(A) W = 1 (16) 

Example 2 

The matrix A is a finite 2n by 2n matrix, and a = 0, thus the sum ([9]) ranges only partitions with 
even number of non-vashing parts. We put 

Aik = (A 2 )ik ■= Si,2c-i , i < k (17) 

Then 

(A 2 ) h =\ l 13 AeSCP ( c ) (18) 
2 ^ 1 otherwise 

where h = (hi, . . . , h^) is related to A = (Ai, . . . , Ajv) as hi = A, — i + N ,i = 1, . . . , N = 2n. 
Example 3 Given set of additional variables t' = (t' 1; t' 3 , t' 5 , . . . ) where we take 

A nm — (A 3 ) nm :— - e ~ u ™- u " Q( n)m )(|t'), a n — (a 3 )„ := e~ Un Q(n){\t') (19) 

Here, the projective Schur functions Q a are weighted polynomials in the variables t! m , Aegt' rn = m, 
labeled by strict partitions (See [18 for their detailed definition.) 

Remark 2. Let us introduce notation t' x = (1,0,0, .. .). It is known that Qhi^t'tx,) = A*(/i) Yli=i 7T"T w ^ ere 

i<j 

Thus for this choice of t' we obtain 

N 1 

(A 3 ) {h) =A*(h)l[— (21) 

z — 1 

One may compare it with Example 5 where f(n) = n. 
Example 4 

Aim — (Ajnm • ^rH-l, m *^m+l,n* (^^) 

Then 

(A.^4 1 iff A = (Ai, . . . , A 2n ) G FP 
10 otherwise 

where /i = (/ii, . . . , ftjv) is related to A = (Ai, . . . , Xn) as hi = A, — i + N ,i = 1, . . . , N = 2n. 
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Remark 3. For some applications we may need further examples. In Examples 5-7 A depends on a given function 

on the lattice denoted by f. In particular one can choose f(n) = n. Below are examples of matrices A whose 
Pfaffians are well-known (see \91$ and references there). 
Example 5 

. . f(n) — f(m) , 

A nm = {A 5 ) nm := /(n) + /(m) ( M ) 

Then for hi = Xi — i + N ,i = 1, . . . , N, we have 

(A 5 ) {h} =A^(f(h)) (25) 

where 



Example 6 

1 - f(n)f(m) 



A - (A \ — tM Z ^( m ) 07-\ 

Sinm — )nm • — ~ ~TJ — \~~F7 \ \^ ' ) 



Then for hi — X% — % + N a i = 1, . . . , N 3 we have 

(A 6 ) {hy =A^(f(h)) (28) 

where 

Example 7 

A = ( A:\ ■ 

(/(n)+/(m))2 



f(hi)-f(hj) 

i<j<N 

A — ( A \ — /( n ) ~ /( m ) 

Anm = (Arjnm ~ — . (60) 



Then for hi = Aj — i + N , i = 1, . . . , N , we have 

(A 7 ) {h} = A$ (f(h)) (31) 

where 



A^v 7) (f(h)) := ( I] 

\i<j< 



N (f(hi) + /(ft,)) 2 ) \f(hi)+f{hj 



m TTTT-rrm ( 32 ) 



Having these examples we introduce the notation 

SP(t,N;U) :=SM(t,N;U,Ai) = (Ai) m e~ u ^ s A (t) , i = l,...,6 (33) 



ASP 
1(\)<N 



In particular we obtain 

S^(t,N;M,U) := £ e"^ s A (t) (34) 



^ 1} (t,iV;{/) := ^ e-^ SA (t) (35) 

Aep 

«(A)<N 

s£\t,N;U,c) := J2 e ~ UXs ^) ( 36 ) 

Aescp(c) 
e(\)<N 

si 1} (t,N,t';U) := ]T e-^Q a(A) (ft') « A (t) (37) 

Aep 

£(A)<JV 

Sf ) (t,iV = 2n,t/) := er u »8 X (t) (38) 

a^fp 

£(A)<iV 

sl 1} (t,N;U,f) := ]T A« (/(/i)) e^s^t) , ? = 5,6,7 (39) 



ASP 

e(\)<N 



The coefficients C^q,} are defined as 



U {a} :=J2U ai , (40) 



G 



The notation U\ serves for 

U x := U {h} , hi = Xi - i + i(X) (41) 

Proposition 1. Sums (|5|). ([M| - (l39|) are tau functions of the "large" BKP hierarchy introduced in [T5] 
with respect to the time variables t. Sums (|37|) are tau functions of the BKP hierarchy introduced in [5] 
with respect to the time variables t'. 

Sums over pairs of strict partitions. In the Frobenius notations [TB] we write A = (a|/?) = 
(a-L,...,ak\l3i,...,/3k). where a = (at, . . . , a k ), at > ■ ■ ■ > a k > and = (/3t,...,0 k ) may be 
viewed as strict partitions. It is clear that £(a) = £(/3),£(f3) ± 1, and we imply this restriction in sums 
over pairs of strict partitions below. 
Now we consider 

S^(t;U,A,B) := 1+ £ e ^-*-i>-^«.> A a s {aW (t) Bp (42) 

Q./3GDP 

where given infinite skew matrices A and B and given vectors a and b, the factors A a and B a are defined 
in the same way as before. 

k k 

U {a} = ^U ai , U { _f } _ 1} = J2 U -Pi-i ( 43 ) 

i=l i—1 

We introduce the following notation 

sl?(t;U) S (2 Ht;U,A t ,A 3 ) (44) 

where i,j = I, ... ,7 and matrices A4 are taken from the Examples 1-7 above. In particular we obtain 
series 

S[f(t;U) := 5>-^* A (t) (45) 
Aep 

Sg\t;U) := 1+ e u i»- u ^s W) (t) (46) 

a,/3eSCDP 

sg\t;U) := 1+ e u m-Ui«, S{a{0) (t) (47) 

qGSCDP^SFDP 

Sif(t,t';U) := 1+ eU{n ~ U{a) Q & ( a )(¥) S(am (t) (48) 

q^gdp 

S$(t;U) := 1+ J2 e Uin - Uiay *Wfi(t) (49) 

qGFDP,/3gDP 

5^(t;[/) := 1+ E ^>-^«>a (a |^(t) (50) 

a,/3GFDP 

:(2)/, .11. 



Sg\t,t'X;U) := 1+ E e c/ {«- c/ {°}Q a(Q) (it') S ( ct | /J) (t)g^ w (it") (51) 

Q,0GDP 

^(t^t/) := 1+ £ e u i»- u ^Q a{a) ^t')s Hp) (t) (52) 

QGDP,/3eFDP 

sg\t;U,f) := 1+ £ e ^>-^«>A«(/(a)) S(a | W (t)A«(/(^)), i, j = 5, 6, 7 (53) 

QGDP.pleDP 

Each Q a (^t') is known to be a BKP [4], [6] tau function. (This was a nice observation of [10], |11|). 
The fact that only odd subscripts appear in the BKP higher times ti m — 1 is related to the reduction from 
the KP hierarchy. 

Proposition 2. Sums (|42l) . (|44|) are tau functions of the "large" BKP hierarchy introduced in [13] with 
respect to the time variables t. Sums (|48p are tau functions of the BKP hierarchy introduced in [5] 
with respect to the time variables t'. Sums (|5Tj) are tau functions of the two-component BKP hierarchy 
introduced in [B] with respect to the time variables t' and t". 
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Remark 4. Let us remind that for the small BKP hierarchy obtained from, KP we have the following [26\ 

S*= A a Q a (t') (54) 

aSDP 

By specification of the data A we obtain 

e- U «*>Q«(±t')> E e- u ^Q a (^t')Q a (it"), £ e- u ^Q a {\t') (55) 

n6DP aGDP aGDP' 

The sums (|55|) are particular examples (see \26^) of BKP tau functions, as introduced in J4%, defining solutions 
to what was called the small BKP hierarchy in flttf . 
The coupled small BKP yields series 

5 5 (t',t",D) := £ Q a (|t')-Da^Q^(|t") (56) 

a,/3EDP 
«(a)=f(/3) 

The coefficients D a ,a in f56|) are defined as determinants: 

D a>f) =det (D ai>f3j ) (57) 
where D is a given infinite matrix. Taking D nm = e c/m_f/ ' l S( n | m )(t) we reproduce ()51|) . 

2.1 Action of ^DO algebra on sums 

Here we shall discribc certain group transformation properties of sums S& and We shall also 

present some Virasoro invariant sums S^K 

Consider the following operator acting on the space of functions of infinitely many variables t = 

(h,t 2 ,...) 

W ( k 9) := res x k (g(D x ) ■ Z(y, x)) \ y=x (58) 
where Z(x, y) is the vertex operator [B] 

Z(y,x) := (eEn^fe— )*n e E» _ {j £ JL_ (59) 

n=0 " 

and where := xd x + \ is the Euler operator and 5 is a given function of one variable. We assume 
that g{D x ) ■ x n = g(n + \)x n . 

The operators wjf^ act as symmetry transformation generators on tau functions, and this action 
may be embedded into the algebra of infinite matrices with the central extension, see [B] and references 
therein. The matrix which corresponds to wjf^ is as follows 

W\f ] := A k g(D) , (A) nm = 5 n , m ^ , (D) nm = (n + <5„, m (60) 

Remark 5. Operators Wjf^ may be also viewed as the element x k g(D x ) of the algebra of pseudodifferential 
operators (^DO) on the circle with a central extension, see AvoendixWi 

In the Preposition below we shall use the notation (^4)_ to denote the antisymmetric part of a matrix 

A. 

Proposition 3. For any data A = (A, a), B = (B,b) where A, B are (infinte) antisymmetric matrices 
and a,b are (infinite) column vectors, consider the following one-parameter family of 



,e tw * e u b 



A(U, t) := (Ve^WeV^) .e'Oa) , B(U, t) := [je^ e u Be u e^) 

(61) 

where U — diag(U n ) and where we assume that matrices and vectors in the right hand sides of equalities 
do exist as formal series in a parameter t. Then 



tw * -S (1 \t,N,U,A) = sV>(t,N,0,A(U, t)), k>0 (62) 



tw, 



* 9> ■ S (2 \t,N,U,A,B) = S i2) (t,N,0,A{U, t),B(U, t)) , k>0 (63) 

where Wjf^ and Wjf^ are given respectively by (|58[) and (|60p . and the exponential e*^ 9> is considered 
as formal Taylor series in the parameter t. 
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2.2 Pfaffian representations 

For 

t=t( X ( M ))=: [ Xl ]+... + [ XM ] 

we have for any N > M = 1 

oo 

S^(t(xi);N,U,A) = a n e~ u -x^ 

n=0 

and for any N > M = 2 we have 



J m>n>0 



U n -U m ( 'm n _ 



Proposition 4. Jbr M = N we have 

S (1) (t(x (M) );W,M) - X^TT pf ^ 

where for N — 2n even 

S i:j = -Sji := (xi - Xj)S {1 \t{xi, Xj),N, U,A), l<i<j< 2n 
and for N — In — 1 odd 



j(Xi -Xj)SW(t(xi,Xj),N,U,A) if l<i<j<2n 
'~ [SW(t(xi),N,U,A) if l<i<j = 2n 



and where 

A N (x) := Yl ( x * ~ x o) 

0<i<j<N 

We shall omit more spacious formulae for the case M =^ N. 
Remark 6. Let us write down the entries of S to express S\ , i = 0, . . . ,6 



S{ 1 \t(x i ,x j ),N,U) = (x i -x j )- 1 e~ u "- u ™W 

m>n>0 
oo 

S{ 1) (t(x i ),N,U)=Y,e~ Vn x? 

n=Q 

oo 

Sj 1 ' (t(xi, xj), 2n, U) = (xi - XjY 1 ^ e - u ^-v c ~ n (x c ~ n xf — x c j ~"x i ) 

n = 

oo 

Sf'(t(Xi,Xj),N,U) = (Xi-Xj)' 1 ^ e ^'*~ Um Q(n,m)(t')(£™^r - Xf3 

m > n > 
oo 

sP(t(xi),N, U) =^e- c/ "Q (n) (t'K 

7l>0 
OO 

sP(t(x,, % ), AT', E0 = (a* - E e~ Un ~ Un+1 i^x n 3 +1 - xfx^ 1 ) 

n>0 

oo 
n>0 

In particular substituting (fT5|) . (|22|) we obtain 



(t, N,U = 0) = Pf 7; r 



/±n(X) 1 — XiXj 



Then it follows that 



and 



N 



x;»A(t(0) = na-^r 1 n a-^) 

AGP i=l i<j<N 



E *AuA(t(x W )) = [] (l-^)- 1 

AGP i<j<N 

Formulae (|5P1) and (15T1) are known, see Ex-s 4-5 in 1-5 of [T5] . 

It is convenient to re-write these formulae in a way independent of the choice of N: 



(79) 
(80) 
(81) 



Proposition 5. 



E s x {t) = ei^i-C+i:^ fe» 



AGP 



X;»AuA(t) = e*E=-^ + E=- 



(82) 
(83) 



AGP 



Relations (|82p and (|83|l toZZ fee used later in Section to solve certain combinatorial problem. From 

SAtr (t) = (-l)l A l SA (-t) (84) 

we obtain 



]T(-l)l A l SA (t) - eiES-i 



AGP 



and 



E SA (t) = e^S =1 -4-E~ =1 

A£P e „ e „ 



(85) 



(86) 



By the simple re-scaling t m — > z m t m in equations 
of z we obtain 



and equating factors before same powers 



Proposition 6. 



E 

iep 
|x|=r 

E (-i) |A| ^(t) 

E S ^UA(t) 

Aep 

A| = r 

E 



Aep 

|A|=T 



S(r)(t) 
s (r) (t) 
s (r)(t) 
s (r) (t) 



|M=T 



(87) 
(88) 
(89) 
(90) 



where auxilary sets of times t = (ii, <2, . . . ) are specified in the brackets to the right of equalities. 
For instance we get ((57)) from l[52|) using the equality 

oo 

E ^ |A| *A(t) = e£™- z 2 ™" 1 *^-* = £ z T S(T)(i) 

AGP T=0 

where t = (t 1} ^,t 3 , ^,t 5 , . 

Formula (l87l) in case t = (1,0,0,...) has an interpretation in terms of total numbers of standard 
tableaux of weight (1 T ) and numbers of involutive permutations of St, see Ex 12 1.5 of [IS], [73], [ST] (see 
also (f2"9l"]l). 

We get from Proposition 2] 
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Proposition 7. 



^ 1} (t(x( 2 ")),2n,;7) 



1 



A 2 „(x) 



Pf 



E 



-U n — U m (mn _ m n\ 



m>n>0 

Choosing U m =0, m < £ + 2n — 1 and C/ m = +oo, m > L + In — 1 we obtain 



(91) 



E ^(t ( x(-))) = _l_Pf 



A 2 „(x) 



(1 - Xi)(l - Xj)(l - XiXj) 



l~( Xl x 3 ) L+N + ^ 



X L + N _ X L + N * 



X n X i 



Proposition 8. 



where 



Si 1) (t(x< a »>),2n,t0 = Pf^-a^/fo^tf)] 2n 

a 2 „(xj :J 



/(z.tf) = J] e-^-^z™ = Sl\t{x uXj ),U) 

m=0 

Choosing U m = 0, m < L + 2n — 1 arid [7 m = +oo, m > L + 2n — 1 we obtain 



E *AuA(t(x( 2 «>)) = -r-^-rPf 



A 2 „(x) 



(Xj Xi) 



1 (x^ Xj ) 



L+2n-l 



1 



i.j — l,...,2n 



Next, as a corollary of Proposition @] we obtain 
Proposition 9. 

1 



x [c+|] -xf + ^ * 



Aescp(c) 



Xj X6 



i 



pf 



-in *J 



/ i,j=l,...,2n 
' c +5l / „ \[c+A] 



i,j= l,...,2n 

(92) 



(93) 



(94) 



(95) 



^ (-l)£?^-i+2«) SA (t(xC 2 «))) 

i,j = l,...,2n 

(96) 

where [a] is equal to the integer part of a. Notice that in case c = n we have only one term related to 
A = and thus the both sides of identity (|95[) are equal to 1 (compare to Lemma 5.7 in f^Vj. 



2.3 Specializations and Examples 

Links with group characters There is a known relation (see [5]) between the Schur functions and 
the odd orthogonal character so\ of rectangular shape as follows 



53 s x (xi,. . . ,x m ) = (xi . ..x m )2 p so//|p\(xf\...,x m ,x^ 1 ,l) 



(97) 



The odd orthogonal characters so\ {xf ,xj , x* , 1), where XJ is a shorthand notation for X\ , x 1 
and where A is an m-tuple (Ai, A 2 , . . . , A TO ) of integers, or of half-integers, is defined by 



so A (x 1 :1 ,...,x* i ,l) 



det 



see, say, (3.3) in J> ; ). Thus, may be equated to a special character of the orthogonal group. 
There is the similar link between S± and a character of the symplcctic group. 



-(Ai-i+m+i) 



det ^.r - - .;■ 



J ' ~ X 3 



(98) 
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Links between sums and matrix models I. Sums as perturbation series for matrix models. 

This topic will be considered separately in [I] . 



Links between sums and matrix models II. Discrete analogs of matrix models. We know [52"] 
few (basically three) ways to choose parameters t in order to convert series in the Schur function s\(t) to 
discrete analogues of matrix integrals where integrals over eigenvalues are replaced by sums over integers. 
These are 

{Al) t = t oc . := (1,0,0,...) (99) 

(A2) t = t(a,x) := a[x] (100) 

(Bi) t = t(g) := (h( q )Mq)Mq),---) , t m ( g ) ~ --^ (101) 

ml — q m 

1 1 - n am 

(B2) t = t(o;g) := (ti(o; g ),t2(o;g),t 3 (o; g ),...) , t m (a;q) := - 1 J 1 (102) 

iV 

(C) t = t(x JV ) := ^ [ Xi ] (103) 

»=l 

The notations [x] and a[x] are standard in soliton theory and denotes Miwa variables 



N = ( x, —, y?--- ) , a N = ^ay^y, ■ • ■ J (104) 

These specializations of t variables will be refereed as respectively the cases (A),(B) and (C) below. 
See Appendix |E1 for cases (j9"9>([rU2|t . 
We have the following observation 

Proposition 10. Let us choose specializations of t variables according to either or (jlOOl) in sums 
© and put 

u n = u ni t*,i*) = up + ± f^fc + tsi» (^) - £ (^V" ^ ( 105 ) 

( \cn + a J \cn + a J ! \cn + a / 

where C/n '' and a,b,c : d : are arbitrary parameters conditioned by ad — be =/= 0. Then sums © are tau 
functions of the "large" BKP, introduced in [T3] with respect to the time variables t* = (t*, t\ , . . . ). 

Proposition 11. Let us choose specializations of t variables according to either (|101[) or (|102l) in sums 
© and put 

" <«™«> " ^ + 1 (££)"t + - t (£3)^ doe) 

m— 1 v ' m—l v ' 

where £/ 9 l°' ) and a, 6, c, d, are arbitrary parameters. Then sums © are tau functions of the "large" BKP, 
introduced in [T3] with respect to the time variables t* = (t*, tjj, . . . ). 

We have 

W 1 

s A (too) = Ajv(^n 7^7 ' h = (h!,...,h N ) (107) 

«— i 1 

where /ij := A; — i + N. 

In case A = (Ai, . . . , A2„) G SCP (see (0|) then in the variables introduced in ([5]) thanks to ([5]) we 
can write 

n 

SA (t 00 )^(A„(y 2 )) 2 n {c _ mc + y . r y = (yi,..-,i/n) (108) 

where y are related to A as 

2/j = Ai — i + 2n — c , i = 1, . . . , n (109) 
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If A is rewritten in the Frobenius notations, A = . . . , a.k\Px, ■ • - > /3fe), then the last relation may 
be written as 

A t (a)A t (ffl ^ 1 

- ns^T^Ti) n sss ,110) 

We also have (see Remark [2]) 

Q a (it 00 ) = Af ) (a) JJ — , a = (ai,...,a fc ) (111) 

1=1 

where 

A N (h):= [] Ai 3) (a):= J] ^-^> (112) 

0<Kj<N 0<i<j<k J 

(A) First we choose the specialization (|99|) . Then the parametrization (|105l) is available. Putting 
t = we obtain from © 

M 

S^t^iV.CU)^ — £ A h(X) e- u i^"^ A N {h) (113) 



AT! 

/ii,...,/iiv— 



in particular 



M N 

S^it^N^^X)) := — IAjvWI J] Mi(^,t*,t*) (114) 



AH 

hi , . ..JiN—O i—1 



^ 1) (t 00 ,7V = 2n,f/(t*,t*)) := - Yl ( A *(y 2 )) 2 n^(2/<,t*,t*) (115) 

yi,---,Vn= i=i 

M N 

S3 1 ' (too, ^V, too, C/(t*, t*)) := 51 A^Af'^JjM^.t*,?) (116) 



hi,...,hN=0 t=l 

M TV 



where 



Si 1 >(t ot) ,iV = 2n,tf(t*,t*)) AjiWne-^^^^Cfci.f.t*), (117) 

^ /li,... : /l„=0 2=1 

^ i (n,t*,t*) = e- v ^' t *' i *^-(n) J i = l,3 (118) 



A(/.) 4 := J] (hi-hjfdhi-htf-l). (119) 

i<j<N 

(compare with [73] where the same expression as the right-hand side was considered in the context of 
random partitions). Here 

n»,f,t-,^ w+ ± (=3)"«. + <f.(=3) - £ (^P" <120) 

m=l v 7 v 7 m=l x ' 

where a, 6, c, d are arbitrary constants, conditioned by ad — be ^ 0. 

Sums (|1 13[) - (|1 17|) . (I122p - (|125[) and (|126p ~ (|129l) may be viewed as discrete analogues of random matrix 
ensembles (compare with [35]), where eigenvalues of matrices are real non-negative numbers. Then (|117[) 
is a discrete analogue of the symplectic ensemble, (|114l) is a discrete analogue of orthogonal ensemble, 
(|115p is a discrete analogue of ensemble of anti-symmetric Hermitian matrices (see Section 3.4 in }17jV 



and (|116p is the so-called Bures ensemble, which describes random density matrices in quantum chaos 
problems, see [36] for the details. 

From double series (|4^|) over Frobenius coordinates of partitions we have 



A/+1 N+l 



^)(too,Ain:=l + £^ £ £ ^' t,) - V, "' t " ) ^7^ (121) 
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in particular 

M+l N+l 



^M-l + Ep E E e^.O-v(o,0 l^a)^! (122) 

k=l^-> oi,...,a»=0 /Si,...,ft,=0 llij=lK + Pj + 1 ) 



6 33 (too, t^, t^, t j- 1 + ^ (ma 2^ 2^ e n fc f , , . 

fc=l v ' Qi,....a fc =0 p U — .Pk=Q lli.j = l\ lJl i ' Pj ' L ) 

(123) 

(1^, t*) := 1+ (124) 



E (fnumz E E 



41 

M+l N+l 7 4 , x a /n\ 

1 e v(/3,t*)- 2 y(a,t*) A^(a)A 2fc (/3) 



S^W) := 1+ (125) 
M+l JV+1 ; 4 , s i4 /m 

1 e 2V(0,t*)-2V{ a ,t*) i ±kW '±k\P]_ 

h « 2fc ) ! ) 2 a uhu^p,h^ nk=i(*+&+ 1 )(*+&+ 2 ) 2 (<*+&+ 3 ) 

Remark 7. Compare with sums obtained in \2_ 

E aeDP A»nti '^f 1 , (126) 
E. eDP A*(a) 2 n- =1 ^^S 11 ' ( 127 ) 

„ 7 -fa ■ (t*)-ff . +1 (t* ) 

E aeDP , A* (a) 4 nt: e ^TCT 1 , (128) 

E^o fe A»A*(/3) nti ^TT (129) 

where we remind that DP' is the set of all strict partitions {pt\, Q2, • • • , <xn > 0) with the property 014 > cti+i+1, i = 

x*/ s4 rr (q»-Qj) 2 ((q»-Qj) 2 -i) , 1Qr ^ 

AW 1J A (a, + «,)'((«, + a,) 2 -l)- (130) 

Interpretation of series (|122|) - (|125|) as discrete version of ensembles of random matrices stays unclear 
for us. 

(B) In the same way the specialization (jlOip yields discrete analogues of circular ensembles in case 
q lies on the unit circle of the complex plane, q = e^^-^, cf> ^ irn, , n € Z. Here the parameterization 

For instance 

M 

SfCt^iV^t*) := ^ E e " n9 ,t * ,t * ) l A ^ ft )l (131) 

hi.....,hjsf— 

M 

^(tfej.^f.t*) := -— J] e_2F(9 ' t%t * ) A «^)' 7V = 2n ( 132 ) 

^ n >' hi,...,h n =0 



where now 



m— 1 v 7 m— 1 v 7 

Remark 8. For g reaZ t/ie correspondent sums may be identified with the so-called Jackson integrals \3$ 
(C) The specialization (|103|) where put x% — e Vi allows to rewrite ([35]) as 

1 M 

= — — V e v ( h '**)4et(e v * ht ) sgaA N (h) (134) 

hi,...,h]y = l 
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which is a discrete analogue of the following two-matrix integral 

J dU J dRdetiTexp ^Tr ^UYU*R + ^ t* m R m \ j (135) 

where the first integral is the integral over unitary matrices and the second is the integral over real 
symmetric ones, dU and dR denote the correspondent Haar measures. Y is any diagonal matrix (a 
source). The matrices are N by N ones. This integral may be viewed as an analogue of the Kontsevich 
integral. 

Then with the same specialization of higher times t we rewrite (|45p as 



M N 



^ = 1 + E E E (sgnA fe (a) ( det 1 ) sgnA fc (^) ) e ^.t')-v(a,t*) (136) 

Each partial sum related to a given k may be considered as a discrete version of the 3-ple integral 

J dU J dA J dB det A n det B n det (£? + UAU^)~ k exp Tr ^ t* m A m - ^ t* m S m (137) 

where f is an unitary and both A and B are real symmetric matrices of size k x k and where dt7, dA 
and g?£? are related Haar measure. 

New hypergeometric functions. Now we specify factors e~ Ux in series (l35l) - (j38|) in order to get 
certain generalizations of hypergeometric functions and basic hypergeometric functions. 
First, we introduce the following hypergeometric series 



p #(<W(a + n;b + n;t) := £ BHr^T ^ ■ Z^ 1 ' 2 ' 4 ( 138 ) 

AeP(3 lli=l l°i + "JA 



£(A)<JV 



and their g-deformed version 



p $f \| n ;b + n ; 9 ,t) := £ flki^l^ SA ( t ) , 0=1,2,4 (139) 



nU(<? bi+n ;<?)A 

« (A)<JV 



where Pi = P, P 2 = SCP, P 4 = FP. 

Here by a + n and by b + n we denote a set of parameters ("indices") (ai + n, . ..,a p + n) and 
(6i + n, . . . , b r + n) where n and N are integer parameters and where the sum ranges over all partitions 
whose length (i.e. the number of non-vanishing parts) do not exceed N . The notation (a)\ where A 
has n non-vanishing parts serves for 



(o)A:=(o)A 1 (a-l)A a --(o-n + l)A B) (o) = 1 (140) 

l) !a 

T(a) 



where (a) m '■— ^Ffffi^ i s the so-called Pohgammer symbol. Then (q a ;q)\ is the g-deformed version of 



(«)a 

(q a -q) x :=(q a ;q)x 1 ---(q a - n+1 ;q)x n (141) 
defined via the q-deformed Pochhammer's symbols: 

(q a ; q) : = 1, (<f ; q) n : = (1 - q a ) ■ ■ ■ (1 - q^ 1 ) (142) 

Remark 9. /n case i m = — 5^ <=1 ^T 1 ^ e summation range is restricted by the condition £(X) < L because the 
Schur functions vanish on the partitions whose length exceed L. 
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Remark 10. As one may notice that by specialization oft = (t, 0,0, ...) and N = 1 we obtain the generalized 
hypergeometric function '381 



p $ r (a + n.b + M)-!, nI=i ( 6i + n)fc W 



Taking t= (f,^, 2 ^, 



iwe obtain 



■which, say, for bi + n = 1 yields respectively the known generalized and the basic hypergeometric functions ]38j. 

The hypergeometric function (|138[) is obtained as a specific case of (|35p if we choose U variables as 
follows 

q p 

U m = ^2 l0 g( & *)m - ^2 l0 g( a *)™ ( 143 ) 

i=l i=l 

The choice 

r p 

u m = iog(9 fcl ; q) m - lo ^ az ' «)«* ( 144 ) 

i=l i=l 

gives rise to (|138|) . 

Hypergeometric functions II We also introduce 



„*(^(a + n;b + n;t):= £ ffHrT^ 



V4 s AuA (t) (145) 



ASP 
^(A)<N 



^)(a + n;b + n;t):= £ ffi* fa + ")a.i (t) (W6) 

lli=l 0* + »ki 



f (A)<]V 



and their (/-deformed version 



p *^(a + »;b + n; ff ,t) := £ ^lg^M^(t) (147) 



ASP 
< (A)<N 



p f^(a + n;b + n; g) t) := £ feg^ *(t) (148) 



1EP™ 
f (A)<JV 



where 



(a) XiP :=(a) Xl (a-^p)x 2 ---(a-^(n-l)l3)x n , (<z)o,4 = 1 (149) 
Thanks to the formula s>(— t) = (— l)' A 's A tr(t) we have 

(a + n; b + n; g, t) = (a + n; b + n; g, -t) (150) 

-(Iff] - (4 JV) 

'(a + n;b + n;g,t) = ' (a + n; h + n; q, -t) (151) 



Remark 11. Having in mind the known relation (see flSf , ]39^) 

f s 2X (XU)dU = 



I uev(N,Wp) - (In) 
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we obtain 



f p ^ N \ a + n;b + n;XU)dU := £ ^ K + ^ ^ (** ) ^ = 12 4 (lg2) 



'(7eD(2fe,F /3 ) 

and their q-deformed version 



ASP 
«(X)<ls 



(Iat) 



C7gU(2fe,F 3 ) 



J (a + n;b + n;g,X[/)dJ7 := £ n^T^T^T^ ' . = 1.2,4 (153) 

1 li = l 



Thanks to results of Subsection 12.21 we have Pfafhan representation for each of the introduced hyper- 
geometric functions. In particular 

p ^' N \a + n;b + n;t) = — Pf (ay - Xi ) p $^ N \a + n; b + n; z^) (154) 

Aat(x) 



p ^^(a + n;b + n;t) - — — - Pf ( Xj - x,) P $r iVj (a + n; b + n; x iXj ) (155) 

Remark 12. iei its note that the hypergeometric series (|138I) . (1145[I . (I145I) are different from the so-called (case 
C ) hypergeometric function of matrix argument 135$ , \34\[ 



E 

AEP 

e (a)<jv 



.Ml 



and hypergeometric series (|139[) are different from Milne 's hypergeometric series ]33f . \34j 



E 

AEP 
« (A)<N 



which are examples of KP tau functions [291, [301. In these formulas H\ and H\(q) are the hook-product and the 
q-deformed hook product respectively. 



3 Fermionic representation 

We suppose that the reader is familiar with the definition of the Fermi operators and the vacuum 
expectation value, for notations see Appendix |Bj 
One may prove the following relations 



S« (t, N, U, A) = (N\ r(t) T(U) g— (A) |0> 



where 



g (A) = g—(A, a) = ^W^+Ensz a n M 

In particular we have 



^ 1} (t,iV;M, U) 
si 1] (t;N, U) 
S^\t,N= 2n;U,c) 

SP(t,N,t';U) 
s£\t,N = 2n;U) 

S?\t,N;U,f) 



(N 
(N 
(N 

{N 
(N 
(N 
(N 
(N 



r(t 
r(t 
r(t 

r(t 
r(t 
r(t 
r(t 
r(t 



e E_M>m>« Tp™^i+T, m <M V>m4>0 |Q\ 

T(U) e^<« | ) 

T(17) e^mez^m-l |q) 

T(C7) e E '»>" JW+zw ^^+^-ez V>m^o | Q ^ 
T(C7) e ^->- /-?W/"n) ^n+Z mSZ </-™<Ao | Q ^ 
T(f/)e E " I >" | ) 



(156) 
(157) 

(158) 
(159) 
(160) 
00 |0) (161) 
(162) 
(163) 
(164) 
(165) 
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where 

r(t) = e ,/ W, f(t) = e J W (166) 

OO OO 

J(*) ■= "^n, J(i) ■= J -™*»' J ™ := V'mV'm+n (167) 

n— 1 n— 1 m£S 

and 

T(E0 := exp ^ t/^Vl - ( 168 ) 

i<0 i>0 

where the fermionic operators are defined as in [BJ, see Appendix [B] 

Remark 13. Pfaffian representation presented in Subsection Ig.gl considered above may be obtained from the 
Wick's rule. 

Remark 14. One can write 

g~~(Al) := e E m >„ *™*n+E m£ z *m^0 _ VK^H^M + ^oJi^rjAc (169) 
g (A2) •= e ^m<n (-l)">2c-ra+l*m _ £ # I _2c_2 <(>(l)+(-l)l!l (170) 

g—(A 4 ) ■- e ^m,eLM m -x = e i iK*- y m*)<i* ( 171 ) 

The corollary of the right hand side expressions is the fact that sums (|35[) . (|36[) and (|38|) may fee re-written as 
certain multiply integrals N-ply integrals for S 1 ^ 1 ' , S^, and N-ply integrals for ), this will considered in 
details in J^. Now, we shall mention a general remark. 

Imagine that a sum @ we can present A nm as moments, or, the same we can solve the following inverse 
moment problem: given A nm = — A mn , m,n > to find such an integration domain D and an antisymmetric 
measure dA(x,y) = —dA(y, x) such that 



A nm = x n y m dA(x,y), n,m>0 (172) 

J D 

Also 

a n — / x n da(x) (173) 
J j 

If we have (|172|1 then in case N — 2n we can write N-ply integral 



St 1) (t,2n,U,A) = e-X?= 1 v* j ( ]~[ , - ; 1 . A,„ (,■ j ) L'f U ,-, ij (.1.71) 
where £ r (t,x) is the following ^DO operator 

OO 

£r(t, as) = J2 *m (xr(D)) m , D = xd x (175) 

m— 1 

and r is related to U as follows 

r(n) = e f/ "- c/ ™+ 1 (176) 
The case U = causes £ r (t,x) = X^m=i ^^"X m and we obtain more familier expression 

2n 

E~ =1 t„» r A2n(a;) pf [(M( ^ Xj)] (177) 



S 1 ' 1 ' (t, 2n, [7 = 0,A)= f TT, 



/n case AT = 2n + 1 u>e ftaue more involved expressions which will be written down in a more detailed version. 
In case the solution of the inverse problem is not unique we have a set of different integral representations for 
the sum @. 

Other representations 

For certain sums like (|161[) we present a different representation as follows 

SP(t,N,t';U) := <JV| T(t) T B (t') T(U) e^->™ *«.(»)*-(»)^*n+E m6 * «c-(m))ifc»*o | ) (178) 



18 



Fermionic representation for . For of (|12")) we have a similar relations 

(t; U, A, B) = (0| T(t) T{U) .9+ (B)g- (A) |0) (179) 

where T(U) may found in (|168[) 

g~ (1) = ff" (A a) = £".™>° ^™^ m +£„ >0 «»tf»*o ( 180 ) 

g+(B) = g+(B, b) = S»,»>o B»«(-i)" +ra *I w _ 1 *i m _ 1 +E„>o (-D^-M-L^ (181) 

For those who are familiar with [T3] these fermionic relations yields a direct proof of Propositions Q] 
and [2] In the next section we will consider it in more details. 

4 Small and Large BKP tau functions 

4.1 A class of "small" BKP (sBKP) tau functions 

We start with the "small" BKP case because it is more simple and illustrative. Concerning this case 
see [26] , 

A general tau function of the sBKP hierarchy may be written as 

T SBKP (t') = (0| r S (t') e £„,m« *»m*»*m |Q) (182) 

and a tau function of 2-sBKP hierarchy as 

T sBKP (t', i') = (o| r B (t') e £».™« A nmMm Tg ^ (183) 

The sBKP hierarchy is obtained from the KP hierarchy by a reduction. In the sBKP reduction, even 
times are set equal to zero and we shall mark it by " ' " : t' = (t^, 0, t' 3 , 0, t' 5 , . . . ), then 

r B (t') = exp J2 r B (tO=exp J2 H -Jn (184) 

n>l, odd 7t>l, odd 

where 

H B = 

2 

The sets of parameters t' and t' are called sBKP higher times. 

"Easy" tau functions We shall consider a simple case, where all terms in the sum of the exponent 
(|182p commute, namely, tau functions 

r sBKP (t',U,A) = (0|r B (t')T B ([/)e£">™>o A„ m 0„0 m +E„ >o «»*»*oVS|o) (186) 

where 

T B (U) =exp (-l)" +1 C/„0„0-„ ) (187) 

We want to single out the [/-dependence though it may be included into the redefinition of A as A nm — > 

P -U m -U„ a 

We shall refer tau functions (|186[) as easy sBKP tau functions. 
We have 

Proposition 12. 

r sBKP (t / ,U,A)= e~ u ^A a Q a (t') (188) 

qGDP 

where the sum ranges over all strict partitions a — {a\, . . . , Ofc), k = 0,1,2, ... , where 

k 



Hn= ( 185 ) 



U w} = u «* ( 189 ) 
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and where Ar a } is the Pfaffian of the k x k antisymmetric matrix A defined as follows: 
for even k its entries are 

A nm = A an0lm , n,m=l,...,k; (190) 

for k odd we take 

A~n m =A anam , n,m=l,...,k; = -A k +i, n = a n , n=l,...,k + \ (191) 

It is assumed that for a = ^4{o} = Qo = e~ u w = 1. 

For proof we notice that e^™>° a ^'t>n4>oV2 _ j _|_ J2 n>Q a n <j) n (j)Q\/2 , and take into account that 

T B {U)4> ai ■••4> ah |0) = e~ ff <->0 ai |0) 
Then we obtain (| 188[) thanks to Lemma [2] in Appendix IE.3I 



Example 1. Choosing 

Anm = Q(n,m)(ji), n > 171 > 0, a„ = Q(„)(|t ) 

where Qr n ,m) is the projective Schur function related to a partition (n, m) and t' are parameters we 
obtain 

T sBKP {t\U,A)= e-^->Q a (t')Qa(tO (192) 

qGDP 

which is actually an example of a tau function (11 83[) , see (31) . 
Example 2. Choosing 

A nm = 1, n > m, a n = l 

we obtain 

T sBKP {t' 1 U)= e- u i^Q a {t') (193) 

Q6DP 

The right-hand side of (|193[) appeared in [69] as a generating function for partition functions related 
to oscillating strict partitions. 

Remark 15. Summation range oo > ct\ > ■ ■ • > ctk > 0, k = 0, 1,2 in sums twer aii strict partitions may be 
replaces by sets of strict partitions whose parts may take values in a given set of natural numbers which we can 
write as a strict partition, say, 7 = (71, . . . ,7at), 71 > ••• > 7jv, N may be infinite. This may be obtained by 
equating e~ Un to zero in case n is not equal to any affi. We obtain the following sBKP tau function 

Tf KP (t',U) := Y, e~ Uia} Q«(t') (194) 

a£DP 

-St 

Relation to solitons. Let us note that for t' = (1, 0, 0, 0, ... ) =: too and 

U n = U n (t*,i*)) = uW -logn! - Y " m C+ Y n ~ m ^ n>0 

m=l,3,... m— 1,3,... 

thanks to Lemma [5] in Appendix IE.4I the right-hand side of (|193[) gives rise to a multisoliton 2-sBKP tau 
function where t* and t* play the role of higher times. Indeed, with the help of (|366[) the right-hand 
side of (|193p reads as 

(0| r B (too) T s (C7(t*,t*)) e^">->o + E»>a |0) = 

OO k 

= i+Y J2 Yl cT " 1 3 ( ° rc ~ ar "' t "~ t/ "° )) j| - - aj - 

= (0|r s (t*)e^«>".>o ^WH e E» >0 e- [ '»%(n)^V2 rB (t*) |o) (195) 

where 



(*) := Y 
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Remark 16. On soliton solutions of integrable equations see \4'tf , ° n multisohtons of sBKP hierarchy see f^. 
The constant U„ plays the role of the initial phase which defines the initial location of the soliton marked by n. 
One can "remove" any soliton by 'sending it to infinity', i.e. via e~ Vn — > 0. See Remark \15[ where 7 may be 
interpreted as the soliton momentums in the N -soliton solution to sBKP hierarchy. 



5 The "large" BKP (1BKP) and "large" 2-BKP tau functions 

A general large BKP and large 2-BKP tau functions may be expressed as the following fermionic expec- 
tation value, respectively : 

T N (l,t) = (N + l\T(t)g\l) (196) 

and 

r N (l, t, t) = (N + l\ T(t)gf (t) |/) (197) 
where N is an integer where T(t) is the same as in (|166p and where 

g = exp ^ A nm lpn^m + BnmVrtV'm + ^'nmV'nV'm + X, ^ n ^ n + 6 »V4)0O (198) 

Here A nm = —A mn and B nm = —B mn . It is due to the presence of 4>q and thanks to equations (|333[) . 
(|33ip the number N may be odd as well as even. 

Parameters t = {ti,t-2, ■ ■ ■ ) and t = (i±,t2, ■ ■ ■ ) are called higher times of the 2-BKP hierarchy. 

The large BKP tau function (|196p was introduced in [TJ] , it solves large BKP Hirota equation written 
down in [13j . Hirota equations for 2-1BKP are written down in the Appendix. 

At the present paper 1BKP tau functions (|196[) are mainly used to study multiple sums, while 2-1BKP 
tau functions (jT97j) will be used to study multiple integrals in pQ. 



Easy tau functions. If for g we choose any product of three special g which are 



g = g (A) := exp ^ A nm ip n ip m + ^ a„V«0o , (199) 

\ n.7n>L n>L / 

fJ = 9 + (B) := (exp ]T iW^-l^m-! + E b nM-n-l) (200) 



n,m>L n>L 

where I is the right hand charge. We will also use 



9 = 9 



-(D) := exp £ D nm i, n ^_ m _ x (201) 



n,m>L 



we obtain simply expressions for (I196p and (I197p . This is because of the fact that all fermionic operators 
in the exponents in (|199p . (|200p . (|20ip anticommute. (In this sense we treat fermions as Grasmannian 
variables). 

The integer L indicating the summation range is arbitrary. Next we consider basic examples. 
First, let us consider 

r=(l\T(t)T(U)g-+(D)\l) = det[D hi , hi ]e- u Ws {h} (t) (202) 



then choosing D nm = S(„| m )(t) we obtain that it is equal to 

c t J2 e- u ^ l) sx(t) Sx (i) (203) 

AGP 

where U\(l) and c; are the same as written down below in Proposition ll3l This is the well-known solution 
of TL hierarchy where the set I, t, t plays the role of higher times. These series were called hypergeomctric 
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tau functions in [55] because they keep many properties of ordinary generalized hypergeometric functions 
(where the role of Gauss equation takes the so-called string equation). Various specifications of this tau 
function were widely used in various problems: in analyze of generalized Kontsevich model |40| . in 
2D chromodynamics [48], [49], [86], some c = 1 string theory calculations [46], evaluation of Hurwitz 
numbers [51] . generalized hypergeometric functions [32] (where the general form (|203l) was studied), for 
models of random partitions [73] . for perturbation series in coupling constants for two- matrix and for 
normal matrix models [21] , |32) . for construction of solvable matrix integrals $T\, for some calculus in 
Seiberg-Witten theory [17] , Gromov-Witten theory [S3], [SI], physics of electronic liquid [S5], models of 
random turn motion by M.Fisher [2T], so-called melting crystals problem [35], [5T] 6- vertex model [62] . 
[63] . [64] . [65] . and in many others problems. 

We hope that the relatives of this series which will be presented below will also find wide applications. 

5.1 1BKP tau functions r N (l, t, U, A) 
We will consider 

T N (l,t,U,A):= (N + l\T(t)T(U)g—(A,a)\l), (204) 

where 

g—(A,a):=ei A ™^„+T, n€Z a„,/>„<Ao V2 (205) 

and T([/) is as in (|168p . Tau function (|204j) vanishes if N < 0. Wc remind that we deal with a pair 
A = (A, a) which consists of an infinite antisymmetric matrix A and an infinite vector a. 
We have the following (compare with Proposition [TJ) 

Proposition 13. 

T N {l,t,U,A) = q e~ Uxil) A {h} (l)s x (t) (206) 

A 

t(\)<N 

where 

N 

U X {1) := Uxt-i+N+l (207) 

i=l 

and A{h\ (0 is the Pfaffian of a matrix A defined as the Pfaffian of an antisymmetric In x In matrix A 
as follows: 

A h (l) := Pf[I] (208) 

where for N = 2n even 

Aij = -Aji := A hi +i,hj+h l<«<J<2n (209) 

and for N — In — 1 odd 

a„ = -a„ := K+«<+* " <2io) 

[a ht+ i if 1 < i < j = 2n. 

We set A Q (l) = 1. 

The constant ci is defined by 

! g-C/i-j U t > q 
1 1 = 0. (211) 

e u l +-+u- 1 l K o 

As we see Ah(0) and U\(0) coincides respectively with Ah and U\ defined in Section [5] This proves 
Proposition [T] 

Remark 17. The right-hand side of (|206|l may be also obtained as a IDKP tau function. 

Remark 18. The right-hand side of (|206p may be also obtained as a special limit of a tau function \225\) below. 
This is a case where e~ Ux vanishes if the length of partition X exceeds N, to ensure this we put e~ u - N - 1 = 0, 
or, the same, we put U-n-i = +oc. 
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As examples we have 

~ lBKp - sBKP (l,t,t',U):= £ e~ u ^ Q l+X - fa) s x (t) (212) 



' N 



7 iBKP-.BKP ( i )t ^ U = Q):= J- Q l+X -fa)s X (t) (213) 

A 

£(\)<N 

r N (l,t,t',U):= J2 e-^ ( °sx(t) (214) 

A 

rjv(t) := J! s ^(t) (215) 

A 

where / + A - denotes the strict partition whose i-th part is equal to Aj — i + N + I . Here t' is the set 
of variables denoted (t'^ttj,^, . . . ) and Q a {\^') is the projective Schur functions [IB] related to a strict 
partition a. As we shall show later 1BKP tau functions (|212[) and (|213D are also sBKP tau function 
whose higher times are t'. 

5.2 1BKP tau functions r N (l, t, U, A, B) 

We begin with a rather special (9(2oo + 1) clement as follows 

g = g+g-, (216) 

where 

g~ = e T, n>m >0 ^■n^ m +Y ln > ^-n't>0 ^ g + _ g E „ > m > ( ~ ) " + V>L m - 1 ^ - n - 1 + £ n > ^° * - n - 1 (217) 

which are exponentials of nilpotents and mutually commuting operators in the Fock space. 

Let us recall that the Fock space T admits a decomposition as an orthogonal direct sum of the 
subspaces Fn of states with charge N 

T = ® Ne zJ r N- (218) 

We have 

g o \0) = |fi>, |0) = ^ Ifijv) (219) 

where each vector |f2jv) belongs to the subspace Fm- The result we need is 
Lemma 1. The vector |firj) is the sum of all basis Fock vectors in Fq: 

\n Q ) = J2 I A) (220) 

ASP 

where A runs over the totality P of all partitions A = (Ai, A2, . . . ), and |A) is defined as 

I A) = (-)*+-+"* Vax • • • Mlfit-i ■ ■ ■ |0> (221) 

where we use the Frobenius notation for partitions f _/#) /.• A = (ai, . . . , afc|/3i, . . . , /3fc) = (ct\j3) where 
ai > ■ ■ ■ > a k > 0, ft > • • • > p k > 0, k = 0, 1, 2, . . .EL 

Then, we obtain the following simplest nontrivial DKP tau function (a version of a "vacuum tau 
function" ) : 

Proposition 14. We have the following DKP tau function 

r (t) := (0| T(t) |fi) = ^2 s x (t) = e^~=i >*m+E~ =1 ta«-i ( 2 22) 

ASP 

where s\(t) are Schur functions. 



3 Here and below fc = will be related to A = 0. 



23 



The second equality follows from the well-known formula [5] 

(0|r(t)|A)=s A (t), (223) 

which is an example of the KP tau function. The third equality in (|222[) follows from the Exercise 1-5-4 
in [18] which should be re-written in terms of power sums. 
Corollary. From consideration similar to [96] we obtain 

esE^^^ + E^viJ-™]!} = | A, I) (224) 

Aep 

Similar relations follows from each equation of Propositions [5] and \6\ 
Now, let us introduce 

T N (l,t,U,A,B) := (N + l\T(t)T(U)g + (B)g-(A)\l) (225) 

where g~{A) and g + (B) are given by (I199p . (l200p . Then by a direct evaluation of vacuum expectation 
value we obtain 

S {2 \t;U,A,B) = To (0,t,U,A,B) (226) 
which is the content of Preposition [5J 

5.3 Mixed IBKP and sDKP tau function 



r := (0|r(t, t', t")e E ".™s° ( D ^<P^+KMl m -i) | ) (227) 

and 



r := (JV|r(t, t' ) e^".™>« D ^^\0) (228) 
^ 53 2- £ ^^ {ctiM Q a (it').s {M (t) (229) 



h!>--->h N >0 =e DP 

£(a)<N 

where 

r(t, t', t") := r(t)r B (t')r B (t"), r(t, t') := r(t)r B (t') (230) 

and where D{ a h } is given by ([57)1 . 

The particular cases are (D nm = e~ Un 5 nm ) 

r(t,t',U)= ^ 2-^ e -^ SA (t)Q A -(if), (231) 

A6P 
l(oc)<N 

r(t,t',U)= 53 2-^ SA (t)Q A -(±t'), ( 232 ) 

< (a)<N 

r(t,U)= 51 e_t/As A(t) (233) 



f (o)<JV 



(t):= 53 SA (t) (234) 



Aep 

l(a)<N 



where A is the partition with shifted parts: A^ := Aj — i + N, i = 1, . . . , N, and where 

N 



U x := 53 u **-i (235) 



i=l 
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5.4 Sums. Modifications of the Schur measure and tau functions. Discrete 
analogs of matrix ensembles 

Here we plan to get use of the considered series in partitions 

Je- ff »«j(t), £ 2-^e-^Q A -(it'), £ 2"^ e -^« A (t)Q A - (if), (236) 

AGP AGP ASP 

and also of the series 

2 e-^ SA (t) SA (t), Yl 2- /(A) e-^Q A -(|t')QA-(|t') (237) 

AGP AGP 

earlier studied respectively in [30] and [3"T] . 

Remark 19. Let us notice that sums over partitions are studied in the context of random partitions. Repre- 
sentation theory Random partitions were started by the school of A. Vershik starting late 60-es and are under 
intensive studies nowadays (see series of papers by A. Vershik, S. Kerov, G. Olshansky, A. Okounkov, A. Borodin 
on this topic). The elegant fermionic approach to this subject was developed by A. Okounkov. 

Let us write down most studied probability measures on the sets of partitions. By ideology developed by A. 
Vershik these are parts of representation theory of linear and symmetric groups. 

The Plancheral measure on the set of partitions of n is defined as 

N 1 

n] UjJ^ II fa-htf, n = l A l> hi:=\i-i + N 

i=l ^ ''' l<i<j<N 

where N > i(X) Q. The z-measure is defined as 

n! TT (z — i + l) Ai (z' — i + 1)a ; tt , n2 , , . , , T 

(i^II (j&5 11 fr-th) , n=\X\, h i :=\ i - i + N 

v ' i=l v ' l<i<j<N 

where z and z' are parameters, and {z)k '■= I 'r<t)' 1 * s ^ ne Pochhammer symbol. 

The Schur measure on the set of all partitions was introduced by Okounkov in [73]. The weight of a partition 

A is 

W x (t,i) = s x (t)s x (i) 

where t and t are parameters of the measure. The Schur measure generalizes the (poissonized with a parameter p) 
Plancheral and z, z'-measures which may be basically obtained as evaluations of the Schur measure respectively 
at the points t = t = e~ p/ Hoc and t = e~ p/2 t(z),t = e~ p/2 t(z') in notations (|38u) . J38TJ. The poissonized 
Plancharel measure on the set of all partitions assign the weight 

, N 1 

e HjhW II (^-M 3 . hi:=Xi-i + N (238) 

i=l ^ *'' l<i<j<N 

The similarity of (|238[) to ensembles of random Hermitian matrices was observed and intensively worked out to 
solve combinatorial problems in late 90-es in papers by Okounkov, Borodin, Johansson and Baik, see [73] , [78] , [80] 
Earlier it was used in physics in [S2], [S3], [S3]. Let us note that evaluation of the Schur measure at other points 
(|382|l and (I383|l yields links with different matrix models, see |32| . 

Now turn out to the studied series in partitions. They bring us to consider the following weights on the set of 
partitions 

W x (t,U) = e- u >-s x (t), W x (t,i,U) = e- u >>s x (t) Sx (t), (239) 

if A = (a\/3) 

Ty A (t,t',t',(7) = 2- f(A) e -^ SA (t)Q tt (it')Q /3 (|t') , (240) 
and, in case partitions are restricted to have at most N parts: 



W x (t,t',U) = 2—^ e-^a A (t)g A - (§*')> (241) 



7 As one can check II238H does not depend on choice of N if N no less than the partition length £(X). 
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where A is a strict partition denned by A = (Ai , . . . , Ajv) as A := (Ai — 1 + N, A2 — 2 + N, . . . , Ajv). On the set 
of strict partitions a £ DP 

W a (t',U) = 2-^e- u i«iQ ct (±t'), W a (t',i',U) = 2' i(a) e- u ^Q a {\t')Q a {\l') (242) 

At last 

Wx(t,U) = e- u ^sxux(t) (243) 

restricted on the set of partitions A U A = (Ai, Ai, A2, A2, . . . ) 6 P 2 may be viewed as the weight of A 6 p. 

Here e~ Ux (or e~ u l"}) plays the role of additional Gibbs-Boltsmann weight assigned to each configuration A 
(or a) induced by external sources. Let us note that in [/-dependence of partition functions Z = Z(U) shows 
non-analytic behavior (phase transitions, compare to |82| . [83) . |86] . |69| . or just to, say, i-Fb(a) = (1 — x)~ a 
behavior which is the simplest example of the series under consideration). Now, WA(t,t,0) is the Schur measure 
studied by Okounkov in [73] while W a (t' , t', 0) is the shifted Schur measure introduced and studied by Tracy and 
Widom in [57] . 

Then all normalization functions (in physics: "partition functions") Z — Wx of these ensembles of 
random partitions are tau functions (|236|) and ([237]). Provided the weights are non-negative the probability of a 
configuration A is 

Wx 

Px = ^ 

where for each ensemble Z is a tau function. 

It is interesting, that if we put all t, t, t', t' to be equal to (1, 0, 0, ... ) and deform U via deformation parameters 
t* according to (|244[1 below we obtain that the partition function Z = Z(f) is again a tau function, in this case, it 
is a tau function with respect to the deformation parameters t*. Moreover, it will be related to discrete versions 
of ensembles of random matrices (where in context of considerations in physics parameters t* are commonly 
called coupling constants). 

The last remark is the following. It is natural to consider bi-measure on pairs of partitions taking general 
tau functions as weight functions. Such models should possess good properties. Tau functions of 2-KP , 2-1BKP 
(©), 2-sBKP (igOOj) and 1BKP coupled to sDKP provide respectively P x P, DP x DP, DP x P with bi-measures 
according to formulae ([3]), ([2]), (|398|) . (|229[) . Objects of these models like correlation functions should be expressed 
in terms of Baker-Akhiezer functions. 

Now we shall write down normalization functions Z for ensembles of random partitions with measures 
W x according to P^l)- (l^]) . Parameters t,t,t',t' will be equated to one of ||55D)> , ([551 ]) . ([552 ]) .([555 ]) 
yielding different specifications of measures (|239l) - (|242p . 

(i) Ensemble with the measure W,x(t, U(t*)). Discrete (3 = 1 ensembles. The partition function 
Z of this ensemble is the tau function (|214[) . 

First we put t = ePtoo taking 

00 

t^ = (1, 0, 0, . . . ), U n = U n (t*) = - J2 ( 244 ) 



Then due to (|384[) the weight function is 



n<m<N 

the normalization function Z — Z^ 1 ' is equal to 



(o) 



zU(t*,p)=T N (t 00 ,U(t*)) = ^ r E II K-h m \ n^'°' h , (246) 

h 1 ,...,h N =0 n<m<N j=l 



where the factor TV! appears when we spread the summation over the cone hi > ■ ■ ■ > > to the 
independent summation over each of hj — 0, 1, 2, . . . at the same time changing Ilr^iT^jvC 1 ™ — ^m) to 
rin<m<Ar ~ h"m\- cn is an unrelated constant. The parameter of Poissonization p can be identified 
with t\. 



2(i 



If we take the same U as in (12441) but take t as in Q38ip we obtain an analog of the Poissonized z- 
measure (with z — a): 

W[ 2) (t*,z) = c N (a) [J (h. 

n<m<N 

with the normalization function 

z^(t*)^r N (t(a),u(n) = c -^ 1 E II K-^U^^hr (248) 

h 1 ,....h N =0 n<m<N j=l J ' 

where cjv(a) is given by (|323l) . 

Remark 20. Boire expressions (|246[l and (I246[l which are BKP tau functions evaluated at points t = too and 
t = t(a) are a/so RftTP toti functions with respect to new parameters t* introduced in (I244|l . See a/so Secfoon 
about interlinks between different BKP tau functions. At the same time these tau functions where t* are higher 
times are examples of the f3 = 1 ensemble where the measure is proportional to a sum of delta functions and 
may be treated as a discrete version of ensemble of random orthogonal matrices with positive eigenvalues, where 
the measure deformation parameters are t*. 



N 

n 



i-iv e T,7=i - ! ' • 



(o) 



(247) 



Next, instead of (I244|) we put 



11 f Ct \ TTL OO 

*" = r^T> U n = U n (t*;q) = uW- £ q nm t* m (249) 



771 — — OO 



By f|383p we obtain, a different model with the weight 

N 

W?\t*; q ,a) = c N (a, q ) J] " fl*") II e^— (250) 

with the following normalization function 

^ 3 )(t*;a,g)=r iV (t(a,g),C/(t*)) = 

= c^) g JT \^-g^\ (ga ; g) ^- +1 eS"^- p ff (251) 

?ii,...,fcjv=0 n<m<N j=l WiDhj 

which, for g e S 1 , may be considered as a discrete analogue of the circulate (3 = 1 ensemble ( [T7]). 
cjy(a,q) is defined in Q323p . Let us choose the limit g a — > (or, the same, t is chosen by (|382[1 ) then 

(q a ;q)h i -n+i -> 1. 

Remark 21 . Expression (|251[) is f/ie _B_R"P tau function evaluated at point t = t(a, g) is a£so a 2-BKP tau function 
(|197|l ratft respect to new parameters {t^, n > 0} and {t^, n < 0} introduced in (|249|l and at the same time for 
q £ S 1 is an example of the fi = 1 circular ensemble \17$ where the measure is singular and is equal to a weighted 
sum of delta functions where the weight parameters depends on t* . 



(ii) Ensemble with the measure W\(t, i, U(t*)) of (I239j) . Discrete /3 = 2 ensembles. The 

partition function Z of this ensemble is the KP tau function (|203|) . This case was considered in the 
paper |32) . here we only write down the most general case where t = t(a, q), t = t(a', q) (which may be 
treated as a q- version of z, z'-measure, with z = a, z' = a') and where U are chosen by ()249[1 : 

Wj»{t*;a,a',q) = 

, \ i i \ TT i h n h m \i TT (g a ;g)fe,-n+i(g a ';g)/ 1 ,-»+i E~=-oo t* n i nhi , OKO \ 
c N {a,q)c N {a ,q) [[ (q - q m ) [[ — e " h i (252) 

n<m<N j=l ((QiQjhj) 

with the following normalization function 

Z^(t*;a,a / ,q) = 
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c N (a,q)c N (a',q) ^ -p-r , hn h ,2 1 r (q a ; q) hi - n+1 (q a ' ; q) ht -n+i £~=-oo t*„Q nhi -U^ 

TV' 2^ 11 W s ; 11 a \ \ 2 6 3 

h!,...,h N =0 n<m<N j=l ) 

(253) 

The last expression which is the 2-KP (TL) tau function (|203[) at the same time is 2-KP tau function 
with respect to the variables {i* , n > 0} and {t* , n < 0}. For g g S 1 it may be identified with a discrete 
version of the one-matrix model of unitary matrices |41] . 

(iii) Ensemble ([Mij) : W\(t, t', t', [/(t*)) where t = t' = t' = t^. 
This ensemble is defined on the set of all partitions A = (a\f3) € P, 

where summation ranges over all pairs of strict partitions a and /3 such that £(a) — £(/3). 

One can convert the summations over the cones a± > ■ ■ ■ > > 0, fi\ > ■ ■ ■ > /3k > to the 
summation over independent numbers on and /3j as follows 



1 j fc! ^,.^=0^.^=0 l\ WW* l\ («< + «i)(A + ft) M «. + A + 1 

(256) 

(iv) Ensemble (J2HJ): W A (t, t', U(t*)) where for J7 and t = t' = see ([233]) , It is defined on 
partitions A whose length do not exceed a given number N. Below hi = Xi — i + N. 



hi + hi V /i,! 

l<i<j<N 1 i=l v 1 



z <«»(t.,=™(t, K .4,= £ n (^n(iv-' 



e 



hi>— >h.jv i<3 "' ' 3 *=1 
hi,- ,/i«=0 i<j J i=l v 7 

Actually this choice of ^4 in (|206l) is exactly related to the Example (|212p which is the tau function of 
mixed lBKP-sDKP tau function. Indeed if t' = := (1, 0,0,...) then it is known that 



A,; Xj 

A,;! Xi + Xj 

i—i %<j j 



see p93p in the Appendix. 

In the spirit of discrete-continuous duality one may expected that there is a continuous counterpart 
to (|258[) . Indeed BKP-sDKP tau function may be chosen as follows 

T^ Kp -' DKP (t, t',U) = (N + /, 0| r(t)r'(t') e f r(t)r'(t') \l, 0) (259) 

= [...[ ^"^ 2 TT^,i,t,t,t',t') (260) 

J J Xi + Xj "~ 

where I, is the Dirac sea level of 1BKP vacuum vector ,t and t are higher times of the coupled 1BKP, 
and t',t' are higher times for coupled sBKP. These "times" play the role of deformation parameters: 

i 00 

d{l(x : t, t, t , t ) == — -=. exp ^ {x t n -\- X ^2n — 1 ^ ^ ^2n— l) 

v 2 „=1 

Similar expressions appear in the study of random density matrices and Bures densities |3.6i • 
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(v) Ensembles (|242|) : W\(t, t' , t', U(t*)) where t = t' = t' = t^. These ensembles are defined on 
the set of all strict partitions a — (a%, a.^ . . . ) G DP. 

Ensembles related to (I242D and their continues versions were considered in [55] as examples of sBKP 
tau functions. Ensembles W a (t',U) where t' = and W a (t',i',U) with t' = t' = too are written 
respectively as 

n 

W (7)( t *) = TT £___££ JT J_ e -^>(*') (261) 

i<j J z— 1 

oo ^ oo n 

*"<*•>- i+Ea e n^ni'- 1 "-"-' c*> 

n— 1 ai,-",oi n =0 i<j J i— 1 

where 

<(«) 

^ i }(t*) = E f7 ^( t *)' ^n(t*):=^ 0) + E n '**> UW = U^ (263) 

i=l i=l,3,... 

and ensembles related to the shifted Schur measure introduced in [57] 

^Hn = n(^) 2 f[(^) 2 e-^ (264) 



. Oii + Otj I xl - \Oii 



2 n / -. \ 2 



^ (8) (t*) = i + e i e n ( ^ ) n ( h ) ^ >(t,) ( 265 ) 



n! ^ V a » + 

The last expression may be identified with a KdV tau function where t* are KdV higher times (compare 
with [85]). Formula f|262[) is a sBKP tau function with respect to variables t* introduced by f|263f) . 

(vi) Ensemble ([2191) : VF A (t,t, f7(t*)) where for [/ and for t = see (gUty , and where 

1 * 

^m(^) = / . X i 
i=l 

The weight 

W^(t*,x) = W x (t 00 ,i(x),U(t*)) (266) 

is defined on partitions A whose length do not exceed a given number N. Below hi = Xi — i + N . 
Then the normalization function is 

oo oo N ^ 

Z^(t*,x)=Y,-- - E (h-hj)l[—x^e- u ^ (267) 

h 1= h N =0 i=l *' 

which may be considered as a discrete version of the so-called one-matrix model with a source: 

N 



f e E~ =1 JVH-C + rrAH dH= f An{z) TT e E «rC+E£x 

J J** fJi 



A, 



where Zj and A^ are respectively eigenvalues of the Hermitian matrix H and the normal matrix A. 
(vii) Ensemble fUg) : TV A (t,t, ?7(t*)) where for J7 and for t = t(oo,g), see ([MU]) . and where 

1 N 

m * — ' 

i=l 

The weight 

W^(t*,x) = W x (t(q),i(x),U(t*)) (268) 
is defined on partitions A whose length do not exceed a given number N. Below hi = Xi — i + N. 
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Then the normalization function is 

oo oo N 1 

zW(t*, X ) = E - E -«*') IT t^-sN -1 ^ < 269 ) 

hl =0 /i„=0 i=l 

which in case g € 5 1 may be considered as a discrete version of models of a random unitary matrices 
with a source: 

S ES— Tr^C + TrAa d[/= I An{z) JJ eEmez e»-iC+EiLx 

J ■ , z i 

2—1 

where e Zi and Aj are respectively eigenvalues of the unitary matrix {7 and the source A. 
(viii) Ensemble ([MTjl : W A (t,t', f/(t*)) where t' = too and where 

1 * 
m * — ' 

i=l 

The weight 

WW(t*,x)=W x {t(x),t 00 ,U(t*)) (270) 

is defined on partitions A whose length do not exceed a given number N. Below hi = A 4 — i + N. 
The normalization function is 

oo , h N '\ 

*">(f,*)= e ^n^^ ( 271 ) 

which is a discrete version version of the following 2-sBKP tau function [26] . 

N 

Z (W) 



7r« z» + *j M 



(ix) Ensembles (|243|) . Discrete /3 = 4 ensemble: Now take W^t, U{t*)) of (|243|) where we recall 
A = (Ai, Ai, . . . , A n , A„) and t is either too or t{q). In the first case, t = too, we have 

2) (t*) = A n (x) JT (272) 

where the set 

Xi := Xi-2i + 2n, i = l,...,n (273) 

A„(x) 4 := [| (a* - Xj f {{ Xi - x 3 f - 1) (274) 

l<i<j<n 

The grand partition function is a discrete version of grand partition function of the symplectic ensemble 

(t*)-C/ aj _i(t*) 



^ n! ^— ' J--*- afi!(ij - 1)! 

ri=0 ii,...,i„ i=l 

In the second case, t = t(q), we have 



" p-t^j (t*)-[^._l(t*) 

^ (t ' )=i - lfl 'n MStr (276) 

where the set of {xi} is the same as before while 

A„(g*) 4 := q- 1 11 - g*') V 4-1 - q*i)(q x < - q^ 1 ) (277) 

l<i<j<7l 

The grand partition function is the grand partition function for a the following ft = 4 ensemble 

00 i 2 p-tMt*)-t^-i(t*) 
Z (i2) = y __ y A n (q x ) 4 TT ' ? (278) 
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6 Interlinks between tau functions of BKP 

Here we will show that for be one of (|M]) - ([TU2t (see also Appendix |R4|) . and for U = U(t*) be 
specified later there exist relations between different 1BKP tau functions as follows: 

(l'|r(tW)T(t*) fl |0 = (/' - l\T(t* + ) g* r(t* ) |0) (279) 

where g* is constructed in terms of a given g according to a choice of t^, and where t*. is the collection 
{t* rn , m > 0} while t* is the collection {t* m , m < 0}. 

Here, in case tW is chosen either via (I38Q|) or via (|381[) . the higher times t* are related to variables 
U by (I105p . namely 



, * i-n-i * ( cm + b\ m „ » / an + b\ v—v { an + b 

U n = U n {t*,t*) = UW + Y — —: C + *oln 3 J - y — ^3 

v ' ' n ^ \cn + dj m u \cn + dj ^ \cn + d 

where a, b, c, d are parameters chosen arbitrary in a way we have no singular terns in the the above sum. 
In particular, if we do not want to have dependence on t* parameters the [/-dependence may be chosen 
as 

U n = U n (t* , c) = - £ ™ m C (280) 

m>0 

as it was done in, say, [5B], or in [57], for the different case, for a TL tau function, while 
In case (jTUTl) . or (ITUSj) by 

v.-u^a.M-um-X (£±|)"c-E (281> 

m>0 v y 7 m<0 v H 7 

In particular case we may take 

U n = U n (t* ; q) = - ]T g nm t* m (282) 

m>0 

as in jS2 or [25] . 

Now we equate a 1BKP tau functions depending on parameters U and thus depending on the param- 
eters t* to a certain multisoliton 2-1BKP tau functions where t* play the role of higher times. Thus we 
present a sort of duality between the U and t variable^. This link may be compared with [55] (see also 
references there). 

Actually this section explains "dualities" between discrete and continuous expressions for tau func- 
tions, see subsection 15.41 Such dualities are also considered in [22] . 

Let us choose (|244|) . Now, in Frobenius notation A = (pi, . . . ,pk\ — q\ — 1, . . ., — qu — 1), we have 

-a(p,Q,o,...) = n gf^fc^ ffl r#TTfrf s ) ^ 

„<m<k ( Pn ~ «»»)(Pm - *») Vn=l T( - Pn + ^(^n) J 

(where numbers p n and q n are related to the Frobenius coordinates as p n = a n and q n = —(3 n — 1 El). 
Written in this form the Schur function may be interpreted as a familiar formula for shift of solitons due 
to interaction where factor in bracket is irrelevant and may be included to the choice of initial position 
of solitons. Therefore 

r(t,[/(t*)) = ^ e -^( t *) SA (t) = 

AGP 

V V TT (P"-P")(g"-g") I A (p» - Qn)- 1 \ E ~ 1 t;( P y-,y) (2 o 4) 



8 Here we convert an observation of |27| to the case of 1BKP hierarchy. 
9 or, more gener 
factor in (2H5) [27) 



9 or, more generally one can take p„ = , <y n = with arbitrary a,b,c,d which keeps linear fraction 
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where {p, q}k means summation over sets of integers p\ > ■ ■ ■ > pk > > qk > ■ ■ ■ > qi which may be 
interpreted as a 1DKP multisoliton tau function, r*(t*), where higher times, t*, are related to U = U(t*) 
of original 1BKP tau function. The fermionic expression for this multisoliton tau function is 

r*(t*) := (0| e J ^ g~(A, p) g+(B, q) |0) = r(t, U) (285) 

where operators g~(A, p),g + (B, q) coincide respectively with g~ (A) , g + (B) if in the last group we replace 
each fermionic Fourier mode Vru respectively by tp(p n ),ip^ {q n )- 



7 Partition functions for certain random processes 

A number of random processes may be described via a sort of partition function Z as follows. We have 
an un-normalized weight W^ V (T) for a transition from a state (J, to a state A during a time interval T. 
To get the probability for this transition, ^-^(T), we divide this weight by the normalization function 
(T) which is equal to the sum of weights to get any state during elapse of time starting from the state 
/z, then, 

jv+a(t) = ^^a(t), z„(t) = hV+aCt) 

where sum runs over all possible states A which may be achieved from initial state /x during a time 
interval T. The normalization function Z provides the condition that the sum of probabilities is equal to 
unity. 

In case the state of a system may be identified with basis Fock vectors (n\ and the transition weight 
W^->a(T) may be written as a matrix element of some given operator o(T) which acts in Fock space: 

WV+a(T) := (MKT)|A) 

then, to get the partition function Z^(T) we need to evaluate 

Z M (?) = (/i|o(T)|fio> 

Examples include random turn vicious walkers model |70j (model (B) in section 4) which may be 
treated as a modification of exclusion processes. Versions of this model were considered in [ST], [75], |80j . 
Relations to 

[88] (section 4.1), see also [89]. In these version walkers move in one direction, then 
Recently, we found a mentioning of the fermionic approach to this problem in the paper by Okounkov 
[50] (Section 1.4.4.). Below we present an explicit evaluation of the probability to get a given final position 
of the walkers in the simplest version of this model where we shall use 1DKP tau function (|222l) . 

Example 1. We are interested in creating a Young diagram (YD) A by gluing box by box in a way 
that each intermediate figure is a Young diagram. It means that at each step we can glue a box only to 
a certain number of admissible places on the boundary of a diagram. A consequence of Young diagrams 
may be called the path connecting initial and final Young diagrams. The number of YD along a path 
will be called the length of the path. For simplicity we take an empty YD (YD without nodes) as the 
initial one. Now one can address few questions: (1) what is the number of paths of length T starting 
at empty YD and ending at a given YD A, (this number will be denoted by Wo->.a(T))- (2) what is the 
number of paths starting at empty YD (empty state relates to /i = 0) of a length T (this number will be 
denoted by Zq(T)). If we consider this creation of a YD as a random process describing the gluing of the 
boxes, then the probability po->a(T) to achieve in T steps a given configuration A is defined as the ratio 
of these numbers: 

Po^a(T) = - 

Let the initial state be given by the vacuum vector (0| and o(T) = Jj, T = 0,1,2,... is a discrete 
time. Then the random process describes the process of creating a Young diagram by gluing boxes to a 
boundary of Young diagram in a way that at each step we get a Young diagram. A single action action of 
Ji on A glues one box uniformly to any admissible place of the Young diagram of A (" admissible place" 
means that a gluing here a box we get a figure which is Young diagram again). Jj glues T boxes one by 
one. Now 

W -»a(T) = (0| Jl\X) (286) 
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is an integer which is equal to the number of ways to create a Young diagram of a given shape A by 
gluing boxes one by one in a way that each intermediate figure is a Young diagram. It is clear that the 
weight of the partition is equal to the duration of time: 

T = |A| (287) 

otherwise the transition weight vanishes. 

The partition function (normalization function) is the sum of all these numbers over final states 
(Young diagrams) A: 

z (t) = (o\Jl\n ) 

KP tau function (|2"2"3")) evaluated at t = ixtoo := (h,0, 0, . . . ) generates all {W ->-a(T), T = 0, 1, . . . }: 
r KP (*itoo) = <0|e**|A> = E f <0|J?|A) = E |^a(T) 

T T 

Let us notice that Wo->a(T) = in case T ^ |A| vanish. Taking into account (|223p we obtain 

Wo^A(T)=T!sA(too)S T ,|A| (288) 

where S is the Kronecker symbol. 

1DKP tau function (l2"2~2"j) evaluated at t = iitoo := (h,0,0, ...) generates all {Z(T) := Z (T),T = 

' 1 '---} ; 

^(txt^) = <o|e^^i^> =E |E w °^(t) - E f ^( T ) ( 289 ) 

T>0 ' A T>0 

On the other hand thanks to the right-hand side of (I222p we have r (titoo) = e^*i +tl which in turn 
is equal to J^t°=o ^i s (i)(^) wriere ^2 '■— (1> 5, 0, 0, . . . ) and where s/ n ) denotes the elementary Schur 
function known also as n-th completely symmetric function |18] . Therefore we obtain 

W T |o2n-T 

W=T! S(T) (t 2 ) = E^-^)T ( 29 °) 

Via saddle point method we find that in T — > 00 limit the main contribution in the sum over n is due to 
n sa I — \ y| which yields for large T 

Z (T) = e^ log T+|(iog2-i)-|v , riogT+0(Vf) 
At last we obtain the answer for the probability to achieve a configuration A in T steps: 

Po^a(T) = ^%4^,|A|, := (1,0,0,...), t a := (1,^0,0,0,...) (291) 

S(T)(t 2 ) 2 

As one can see the probability to achieve the state A = (T) is given by 

Po-Kt)(T) 



^o(T) 



One can ask, given T, what is the configuration A which maximizes the number Wo-va(T) and thereby 
the probability po-* a (T)- The answer is known 21 and is given by Kerov-Vershik formula for the so-called 
limit shape of Young diagram, see [71][72 0. 

If we modify our random process and admit both creation and elimination of a box we obtain basically 
the same shape of A however now the weight of this configuration will be less than T [21] : 

Example 2. Starting from the vacuum zero Young diagram, at each time step we either add or 
remove a box at random in a way that a figure we obtain during at each time step is a Young diagram, 



°The reason of this coincidence is the following. The Vershik-Kerov limit shape YD maximizes Plancheral measure on 
partitions, and Plancheral measure is basically (s,\(too)) 2 . 
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see the figure below. This model is equivalent is a model of random turn walk suggested in [70] where 
initial configuration of walkers is the step function. 




1. Random turn walk 2. Random adding/removing a box to a Young diagram 

of particles related to the up/downward hops of particles on Maya 

on a Maya diagram. diagram. At unit time instant either a box has to be added 

at any of vacant places marked by star, or a box marked 

by x has to be removed. 



This model describes hard core particles ("walkers", "hard core" means that two particles can not 
occupy the same site) situated at the sites of ID lattice. The model implies that at each tick of clock one 
chosen at random particle hops either to the left or to the right. In our picture particles are fermions 
which we placed on the vertical lattice. The step function is the vacuum state, (0|, describing the Dirac 
sea where all sites downward to the sea level are occupied. We count particles from the top. Excitations 
may be described by partitions A = (Ai, A2, . . . ) where Ai describes the upward shift of the up-most 
particles (the particle number one) with respect to its original position in the Dirac seas, the shift of the 
particle number i is equal to A^. It is clear that Ai > A2 •• • . Such configuration is denoted by (A|. Each 
upward step of a particle from a configuration A may be described as gluing a box to the Young diagram 
A, while each downward step is described as removing one box from the Young diagram of A. 

The number of paths of a length T which start at the vacuum configuration and end at a given 
configuration A divided by the number of all paths of the length T which start at the vacuum configuration 
defines the transition probability po->a(T). 

Now the number of ways to achieve a given configuration A during a lapse of time T starting from 
the vacuum (step function) configuration is 

Wq^a(T) = (Q|(Ji + J-i) T |A) (292) 

where T is not necessarily equal to |A|. 

The number of ways to achieve any configuration in T time steps starting from the vacuum configu- 
ration is 

Z (T) = <0|(J! + J_!) T |O ) 

A usage of Baker-Campbell-Hausdorff formula e 4Jl+tJ - 1 = e^' 2 e tJ ~ 1 e tJl may be considered as an advan- 
tage of the fermionic approach. After some algebra we obtain [5T] 

Proposition 15. We have 

W ^ x (T) = fi-s x (t oo )6(T,\\\), 6(T,\X\) :=2^ — L- (293) 

I 2 )■ 

where S(T, |A|) replaces the Kronecker symbol 6 t,\\\- I n (|293p T— |A| is an even number otherwise Wo-+\( f) 
vanishes. 
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The following 1DKP tau function generates partition functions 

1 1°) = E S E hwt) = E S z w ( 294 ) 



i e tji+tj- 



T ! Z^, "«-™v-/ T |" 

T>0 A T>0 



On the other hand 



| e *Ji+tJ_i| n j = e5 * (oje^in) = e* +t = ^ < T s (T) (t' 2 ) (295) 



T=0 



where t 2 := (1, 1, 0, 0, ... ). Therefore we get 



Proposition 16. The number of paths of the length T which start at the vacuum configuration is given 
by 

Z (f) = 7ls w (t 2 ), t' 2 := (1,1,0,0,...) (296) 
where S(r)(t 2 ) is the elementary Schur function related to the partition (T): 

[I] 

T7ie number of paths of a length T which start at the vacuum configuration and end at a given configuration 
A divided by the number of all paths of the length T which start at the vacuum configuration is given by 

po^x(f) = ^^S(T,\X\), ti:= (1,0,0,...), t 2 := (1,1,0,0,0,...) (298) 

As one can see 

*WT) = ^y 

For large T one may apply the saddle point method to evaluate the sum (|297[) . The saddle point is 
related to n = § - \<J\+0(1). This yields 

Proposition 17. In large T limit we obtain 

Z {T) = ei log T+ i log I +°(^ (299) 

Po^x(f) = S A(ti)e-^ 2 +°^ (300) 

'One dimensional dimer' target configurations At last let us focus on the following problem. 
Let us evaluate the number of paths of a given length T which end on a configuration related to a fat 
partition A U A : 

A^fp(T) := ^aua(T) (301) 

ASP 

The number TVpp (T) vanishes in case T is odd. Let us recall that each configuration A U A describes a 
configuration of pairs of particles ('one dimensional dimers'). 

Proposition 18. 

, r (2^ T (T- 1)!! iff T is even . . 

N FP (T) = 1 y ' JJ 302 
V ' [0 iff T is odd X ' 

Indeed, on the one hand 

£ (0|e tJl+t ' /_1 |A U A) =E^E ^aua(T) - £ ^fp(T) (303) 
Aep T>0 ' A T>0 



On the other hand (see ((82 



(Q\e tJ i +tJ -i | A U A) = e**" £ (0|e 4 Jl |A U A) = e* 2 = £ (304) 

A6P AGP T=0,2,4,... V2 )• 
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A Appendices 

A.l Pfaffnians. Partitions. Schur functions 

(A) Pfaffians. We need the notion of Pfafhan. If A an anti-symmetric matrix of an odd order its 
determinant vanishes. For even order, say k, the following multilinear form in < j < k 

Pi [ A ] : = S S n ( CT ) At(1),<t(2) Ar(3),a(4) ' ' ' K{k-l),cr(k) (305) 

(7 

where sum runs over all permutation restricted by 

a : a(2i - 1) < a{2i), cr(l) < <r(3) < • • • < a(k - 1), (306) 

coincides with the square root of det^4 and is called the Pfafhan of A, see, for instance [T7]. As one can 

see the Pfafhan contains 1-3-5 (k — 1) =: (k — 1)!! terms. 

The following equality is known as Schur identity 

IT (f^—^ l=A$n(*) (307) 

where 



Xi + Xj y l<i t j<2 n 



A£,(*):= n ^rr 1 ( 308 ) 
i<i<j</t 1 j 

Let us mark that a special case of this relation is obtained if xi n vanishes. In this case we write 

Pf(A)=A5 n _ 1 (ar) (309) 
where A is an antisymmetric 2n x 2n matrix defined by 

if 1 < i < j < 2n 



, = J J (310) 

Hafinans The Hafnian of a symmetric matrix A of even order JV — 2n is defined as 

Hf (A) := A .( 1 ) iCT (2)A CT ( 3 ) iCT (4) • • • A CT (2ri-l) )0 -(2n) (311) 

where sum runs over all permutation restricted by 

cr : a(2i - 1) < a(2i), ct(1) < <r(3) < ■ ■ ■ < a(2n - 1), (312) 

As one can see the this sum contains 1-3-5 (27V — 1) —: (2N — 1)!! terms. 

Remark 22. Let us note that entries on the diagonal of the matrix A does not contribute the sum (|311|l . 
The following equality was found in |91j 

pJ(^i^) )= TT x -^m{(^—\ ) (313) 

+ J^Ktn) l<i<J<2k Xi + X i \\ X * + X iJl<iJ<3n) 

Another proof of this relation was presented in [90] . Let us mark that a special case of this relation is 
obtained if xi n vanishes. In this case we write 

Pf(2?) = A^Or) Hf(C) =: A%(x) (314) 

where B and C are respectively antisymmetric and symmetric 2n x 2n matrices whose relevant entries 
(see Remark |2"2"|) are given by 

U if i<j = 2n, 13 U i/ i<j = 2n, 
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(B) Partitions. Polynomial functions in many variables, like the Schur functions, are parameterized 
by partitions. 

Let us remind that a partition of certain number n is an ordered set of integers A = (Ai, . . . , Aj) where 
Ai > • • • > A; > such that n — X)l=i ^i- Then, n is called the weight of A and commonly denoted by 
|A|), see [TS]. Integers Afc are called parts of the partition A. The number of non-vanishing parts of A is 
called the length of A and will be denoted by ^(A). 

Almost everywhere throughout the paper we will denote partitions by Greek characters. 

Strictly ordered sets a — (a\, . . . ,ai), a\ > ■ ■ ■ > ai > are called the strict partitions, see |18j . 
In this paper we write {a} to denote strictly ordered sets a\ > ■ ■ ■ > ai where numbers a k are not 
necessarily positive. 

We basically use notations adopted in [18] . 



Young diagrams The (Young) diagram of a partition A is defined as the set of points (or nodes) 
(i, j) € Z 2 , such that 1 < j < A.;. Thus, it is a subset of a rectangular array with £(X) rows and Ai 
columns. We denote the diagram of A by the same symbol A. For example, 



(316) 



is the diagram of (3, 3, 1). The weight of this partition is 7, the length is equal to 3. 

The partition whose diagram is obtained by the transposition of the diagram A with respect to the 
main diagonal is called the conjugated partition and denoted by A*. 

In the Frobenius notations (see [TB]) we write A = (ax, . . . , ak\/3\ ■ . . , /3k) or just A = (a\(3). For 
instance the partition (|316p is written as (2,1|2,0). The partition X tr :— (f3\a) is called transposed to 
X=(a\/3). 

Hook polynomials, Pochhammer symbols The product of hook lengths H\ is defined as 

H \ = II h v = A, - i + A} - j + 1 , (317) 

where the product ranges over all nodes of the diagram of the partition A. 
Given number q, the so-called hook polynomial H\(q) is defined as: 

Hx(q) = J! ( X - ^ )' hi = Xi—i + Xj—j + l (318) 
i,je\ 

In what follows, we also need notations: 

fe 

n(X) ;= Y,(i ~ l)Ai , (319) 
j=i 

r(a + m) 

(a a := i a, (a - 1 a 2 • • • (a — « + l)\ k , (a m := — , 320 

T(a) 

{q a -,q)x-.= {q a ]q)xM a - 1 i<l)M---{q a - k+l \q)x k , (<Z a ; q)m ■= (1 - q a ) ■ • • (1 - q^ 1 ) , (321) 

where k = ^(A). We set (a)o = 1 and (q a ; q)o = 1. 
Useful relations are 

k k 

{a) x = cfe(a) JJ(a)h i+ i_ fe , [q a ] q)\ = c k {a, q) J| (q a ; q) hi +l-k (322) 

2 = 1 i = l 

where 

fe k 

Ck(a) = ~[[(a-i) k -\ c k (a,q) = ~[[(l-q a - t ) k -\ h i = X i -i + k (323) 

i=l i=l 
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Schur functions. We now consider a semi- infinite set of variables t = (t\, t%, . . . ). Given partition 
A, the Schur function s\(t) is defined by 



s A (t) = det(/i Ai -i+j(t)) 1 < ii <^ (A) , where ^lz k h k {t) = exp z m f TO , (324) 

fc=0 m=l 

and, for k < 0, we put hk = . The /ifc(t) is called the elementary Schur function. 

There is another definition of the Schur function; it is the following symmetric function in the different 
variables x :— :— (xi, . . . , x n ), where n > £(\): 

det(x Xi ~ j+n ) 

*x(*) = TT7 n—j\ n ; (325) 

det \x 4 L 



for the zero partition one puts Sq( 2 ') = 1- If 



t = t(xW) = (t!^),^^), . . . ), t m (x^) = - 



m 

i=l 

then definitions (|324p and (|325[> are equivalent [18]: 

»A(t(i (n) ))=| A (i (n) ). (326) 

Remark 23. From definition [32$ it follows that s\(t(x^)) = if > n. 

The Schur functions s^(x%, . . . , x n ), where £(X) < n, form a basis in the space of symmetric functions 
in n variables. We use the underline in s x only to distinguish the two definitions. If an n x n matrix 
X has eigenvalues X\, . . . , x n , we may denote s x (xi, . . . , x n ) by s\(X), without underline, since in this 
paper the Schur function with uppercase argument is used only in this sense. 

One of the wonderful results of Kyoto school is the formula 

8 X (t) = (0|r(t)|A) (327) 

which may be obtained by the direct calculation using (|324[) and r(t)^;„r(t) _1 = E™=o'/'»-m' l mW- 

We want to mark out that apart of traditional notation s\ we shall use also notation s^y widely used 
in physical literature, say for instance [S3]. Here hi = Aj — i + N where N is the length of a partition A. 

B Charged and neutral free fermions. 

For the charged fermions we shall use notations and conventions adopted in [5] . In particular for charged 
fermions we have 

WL n -il = W„=0, V-n-i|0)=Vi|0)=0, n> (328) 

^rni'n + tynVL = S nm , n € Z (329) 

Such fermions are used to construct tau functions of KP, Toda lattice (TL) and the large 1DKP hierarchies. 

For "small" BKP hierarchy and "small" BKP hierarchy coupled to 1DKP we need neutral fermions, 
{4>n}nez, see [B], defined by the following property: 

4>n4>m + 4>m4>n = (-l) n 5 n _ m , n£Z (330) 
<Pn1pm + 1prn<Pn = 0, 4>n4>L + V'm^n = ( 331 ) 

It results from (|330[) that (4>o) 2 = |. 

The action of neutral fermions on vacuum states are defined by 

0„|O)=O, (O|0_„ = O, n< 0, (332) 

4>o|0) = -J=|0), (o|0 o = i=(O| (333) 
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Note that the action of 4>o on the vacuum vectors is different from the action defined in [5J . This causes 
a modification in formulation of the Wick's relations, see below. Here we follow [13], where the choice of 
the corresponding Fock space is different from the suggested in [B]. See also the Appendix in the arxiv 
version of [69] for some details. 
For n, m € Z we have 

(O|Vn0 m |O) = (0|^ m |0)=0 

The right vacuum vector of a charge I is defined via 



• - • Vo|0) if/>0 
^••>ii|0) if/<0 



The dual vector is defined as 



(l\ = \?M~*U *'>0 (335) 
U \ (0|^„! • • - -0; if Z < K ' 

Basis Fock vectors. We shall use the following notation 

| A, /) = (-)*+-+* ^ ai+l • • • ' ' ' W-fc-x 10 ( 336 ) 

where A = (ai, . . . , afe|/3i, . . . , /3fc). Dual vectors may be defined via (A',Z'|A,Z) = Sx,\'Sw. The vector 
(|336p may be also written as 

|A, I) = ip hl +i-N ■ ■ ■ ^Ph N +i-N \l - N) (337) 

where A = (Ai, . . . , A at) and 

hi:=Xi-i + N, i=l,...,N (338) 
which are called shifted parts of A. Here N > fi\ . 

Wick's relations. Let each of Wi be a linear combination of Fermi operators: 

Wi = ^ w im^m + W tm^L + Y W im4>m > i = h ■ ■ ■ i ™ 

m£2 m£Z mGZ 

If n is even 

(l\wi ■ ■ ■ w n \l) = Pf [A] (339) 
where A is n by n antisymmetric matrix with entries 

Aij = {l\w l w j \l) , i <j 

Now turn to the case where n is odd. In case zero mode operator 4>o is absent in series for Wi, namely 
each Wio = 0, we have the standard situation where (l\wi ■ ■ ■ w n \l) vanished^. 

C Wick's rule and Pfaffian representations 

From the fermionic expression for tau functions we obtain its Pfaffian representations as a result of 
application of the Wick theorem. 

First we suppose that g may be factorized into the product g — g~ g + such that 

g = g-g + -- (i'\g- = {l'\, 3+10 = 10 (340) 

Let each Wi, i = 1, . . . , 2m is a linear combination 

nsEZ n€Z 



11 We exclude <f>o because we use Kac-de Leur Fock space |13| where (0|</>o|0) = 1 non-vanishes, see Appendix lE.il 
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If n + m = 2k by Wick theorem we have 

• • -w n gw n+ i ■ ■ ■w 2m \l) = Pi[W i:j }\ij = i t ... t 2k (341) 
where W is the 2k x 2k antisymmetric block-structured matrix with entries given as follows 

Wij = (l\wiWig\l), W ip = (l\wigw p \l), W pq = (l\ g w p w q \l) (342) 
where l<i<j<p<q<n. 

D On the central extension in Lie algebra of ^>DO 

The central extension in the algebra of ^DO may be chosen via the choice of nontrivial 2-cocycle in this 
algebra. It is known that there two independent cocycles found in |95j : 

ujAA, B) = res res (A [log x, B]) , u 2 {A,B) = res res (^4 [log d x , B] ) 

x dx x 

In our case the choice of the cocycle should reproduce the Japanese cocycle in the algebra of infinite 
Jacobian matrices [6] and it should be chosen as follows 

u(A, B) = res res (aT^log D x , B] ) , A,B £ ^DO 

x D m 



where the commutator of logD x with a ^DO is defined via relations 

[In A, , f{D x ) \ = 0, [log D x , f(x)] := f) ^ ' f{ * )] °' 

n=l 

Up to coboundary term w = u>i + w-i- 



E Hirota equations 

Fermionic form of Hirota equations was invented in the papers of Kyoto school, see for instance [6] and 
references therein. 
Introduce 

Si := i-) n <t>n ® <t>-n , S 2 :=J2 H n 4>n ® 4>-n (343) 

•n^O n^O 

S := <t>o ® (344) 

S 3 :=Y,ip n ®i>i, Si := ^n® i>n , (345) 

Hirota equations in the fermionic form are: 
For KP hierarchy [B]: 

\g®9,S 3 } = (346) 
For large (the same fermionic) DKP hierarchy [T5] : 

[«?® 5 ,S3 + S 4 ] = (347) 

For large (fermionic) BKP hierarchy [15] : 

® g, S 3 + S 4 + S ] = (348) 
For small DKP hierarchy [4] (see also [13]): 

\ g ®g,S 1 }=0 or\g®g,S 1 ]=0 (349) 
For small BKP hierarchy [3] (see also [15]): 

[ 5 g, Si + S ] = (350) 
For two-component BKP hierarchy in form |13j : 

[g ® .g, Si + S 2 + So] = (351) 
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E.l A remarks on BKP hierarchies |13j and [6 J and related vacuum expec- 
tation values 

Let us note that different vacuum states were used in the constructions of BKP hierarchy in versions [T5] 
and [5J. If we denote the left and right vacuum states used in [5J respectively by '(0| and |0)' then 

<0| = ±= '(0| + '(O|0 O , |0) = -^=|0)' + O |O)' (352) 

Introduce also _ _ 

'{l| = V2'Wo, = V20 O |O)', (353) 

then '(0||0)' = = 1 and instead of (|552l we have 

<i>nW = M 1 )' = 0, '{0\<P-n = '(l\<t>-n = 0, n<0 (354) 

see [5J for details. 

Correspondingly Fock spaces used [12] and [B] are different. From the representational point of view 
this definition is somewhat more convenient, since each Fock module remains irreducible for the algebra 
Boa which is the underlying algebra for KP equations of type B (BKP), see |13) . 

The vacuum states '(0| and |0)' are more familiar objects in physics. In particular any vacuum 
expectation value of an odd number of fermions vanishes, while, for instance, (0|(fo|0) = 4j. 

Let F be a product of even number of fermions. Then it is easy to see that 

(0|F|0) ='(0|F|0)' (355) 

E.2 A remark on formulae containing Q\ functions 

From |10j it is known that 

'<oie*ew* -^io)' = jr; Q( r 'hh n)(§) for N even ' ^ 

10 for N odd, 

V^'(0|^^^---^JO/=(f f f Q ^--^)^ for N odd, 

10 for N even, 

where Q\ (§) , A = (Ai, A2, . . . ) are the projective Schur functions, see [T£]. Thus 

(0\e H ^Mx 2 ■■■4>x N |0) = (-±='<0| + '(0\<t>a\ e H ^Mx 2 ■ • ■ <f>x N (-^W + 0o|O)' 

i '{Q\e H( ~ s UxA^ ■■•fa* W + 'm^ H{s) Mx, ■ ■ ■ faM)'+ 
+ _L '<O|0 oe H «^ A2 ■ • • fa N |0)' + i '(O|e^0 Al A2 • • • faMY 

= l'(0\e Bis) Mx 2 ■ ■ ■ fa N \0Y + ~ '{l\e H ^Mx 2 • ■ • <f>x N \l)' 
+ '(O|0 o e ff(s WA 2 • • • <M0>' + '{l^^UxAx. ■ ■ ■ <t>x*W 

= 2 T Q(Ai,A 2 ,...,Aw) (g) ' 

since the role of ' (0| and '(1| (resp. |0)' and |1)') is interchangeable. 
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E.3 Projective Schur functions (Q-functions) and neutral fermions 

There are two bases of neutral free fermions 

<fc = ^(Vi + (-!)>-,), k = -j=tyi - (-l)V-i), (358) 

where i € Z, each of which generates this subalgebra. 

Using the results for charged free fermions, the anticommutation relations are 

&]+ = 4>j}+ = 4>j]+ = 0, (359) 

and, in particular, 4% = = \. Similarly, the vacuum expectation values of quadratic elements are 
given by 

U-l)%,-j ^<0 

'(0|<^|0)' =' <0|^4-|0)' = Hs jt0 2 = 0, (360) 

[o i > 

and Wick's Theorem is used for arbitrary degree products. 

The neutral free fermion generator is defined by <fr(p) = X^nezP™^"- 

We have (for \p\ ^ \p'\) 
'(0\^p)^p')\0y = l^-l, (361) 
and '(0\4>(p')(l>(p)\0y = -'(0\<t>{p)<j)(j(/)\0y. By Wick's Theorem we get 



2 -n/2JJP^1 TV even 

i<i^ + Pi • (362) 

otherwise 



The connection between the charged and neutral free fermions can be expressed in terms of the 
generators as 

- qiP(pW (-q) + pi>{q)i>* {-p) = Hp)Hq) + 4>(p)4>(q)- (363) 

In the sBKP reduction, even times are set equal to zero and we define t' = (i' 1; 0, t' 3 , 0, t' 5 , . . . ), and 
the hamiltonian 

H B (t')= ]T H B t' n , (364) 

n>l, odd 

where 

ff^Ec-iy+W-i-n. (365) 

For the fermion generating function one has 

0(p)(t') = e HB ^'U(p)e- HB(t,) = e^ B(t 'U(p)e-" B(t,) = e^'Uip)- (366) 

Note also that 

H(t') = H B (t') + H B {t'), [H B (t'),H B (t')} = 0. (367) 
Similar to the KP case, sBKP r-functions are defined by 

r B (t') = (h(t')), (368) 

where h is the Clifford algebra of the neutral free fermions fc. The n-soliton r-function is obtained by 
the choice g = exp(^™ =1 <k<f>(pi)<j>(Qi)) ■ 
The Schur q polynomials are defined by 

exp(2£(p,t')H5> fe (tV- (369) 

fe>0 
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Thus 

^(t / ) = X)%(i t/ )^-fe- (370) 

fc>0 

We have 

'(0|^(t')^(t')|0)' = i % (it')g,(|t') + ^(-l)Wi(|t')*- fc (|t'). (371) 

fc=i 

Since 

1 = exp(2£(p,t'))exp(-2£(p,t')) = ^ q t (t')q^(-t') = ^(-1)^ q^q^t' )p\ (372) 
for all n > we have 

n 

]r(-ir ft (t')<z„- 4 (t') = o. (373) 

i=0 

This is trivial if n is odd and if n = 2m is even then it gives 

m 

< Zm (t') 2 + 2^(-l) fe g m +fc(t')'Zm-fc(t') = 0. (374) 

fc=i 

We can also define 

b 

««,&(t') = ga(t')ft(t') + 2^(-l) fc (Za+fc (t')<Z6- fc (t'). (375) 

jfe=l 

If follows from the orthogonality condition (|373[) that 

g«,&(t') = -*,«(*'). (376) 
and in particular, q a ,a(t') = 0. Comparing (|371|> and (|375[) . it is clear that 

g 0l 6(|t') = 2(0 a (t')^(t')>. (377) 

Now consider A = (Ai, A2, • . • , A2n) where Ai > A2 > ■••A2n-i > ^2n > 0. Note that this is a 
partition with an extra trivial part included if necessary to ensure that the number of parts is even. 
The set of such strict, or distinct part, partitions is denoted DP. For A € DP we define 

Q A (it') = Pf( 9 A 1 ,A J (it')). (378) 

This is the Schur Q-function. By Wick's theorem, we come to 

Lemma 2. 

Q A (±t') = Pf(2(0 Ai (t')^(t'))) = 2"(<K 1 (t>A 2 (t') • • -0A 2 „(t')). 

This wonderful result was obtained in [10] (see also |11) where it was independently found that Qa is 
an example of sBKP tau function). In particular, from this Lemma and (I358P . (|327p it was obtained [IU] 

2-^Q A (it') = v /^0, (379) 

where S\(t) is the Schur function, and the partition A 6 P is the double (see l,I,Ex9(a) of [T5]) of the 
strict partition A. £(X) is the length of A (the number of non- vanishing parts of A). 
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E.4 Schur functions and projective Schur functions evaluated at special val- 
ues of its argument. 

We introduce the following notations: 

too = (1,0,0,0,...) , (380) 
t(a,l) = (|,|,|,...) , (381) 

t(oo,<z) = (ii(oo,g),i 2 (oo,g),...), t m (oo,q) = _ , m=l,2,..., (382) 

1 - (q a ) m 

t(a,q) = (ti (a, q),t 2 (a, q),...) , t m (a,q) = _ , m=l,2,... (383) 

For various purposes these choices of times were used in [33] , [3U] , [2Z] j EH] , [SZ] j [SS] , |H2] • Note that 
t(a, (?) tends to t(oo, g) (resp. t(a, 1)) as a — > oo (resp. q — > 1). As for too, if / satisfies f(cti,c 2 t2, c 3 t3, . . .) = 
c d f(ti,t2,t3, . . . ) for some d £ Z, we have h d f(t(oo,q)) — > /(too) as S := lng — » 0. 

Lemma 3. For a partition A = (Ai, A2, . . . ), let hi := n + Aj — i (1 < i < n), where n > ^(A)). Then 

SA(too) = = 1^-7 , (384) 

SA (t(a, i)) = — = g r(a -z + i) ' (385) 

a «(A) A(VM 
SA (t(oo,g)) = j—r = g L , (386) 

s x (t(a, q)) = g " ( ^ (g . a; . g)A = Cn(a, q) f[ (? a ; q) hi -n+l , (387) 

n n n 

i<j t<j i=l 

where for # A , /f A (g), n(A), (a) A and (q a ;q)x see respectively (j3T7)) . (j5T5)t . J3I2J, (JSSjJ and ([52T]) . Note 
that those quantities (|384l) - (|387[) are independent of the choice of n > £(X). 

Lemma 4. . Let (a\f3) = (a±, . . . , (Xk\Pl> ■ ■ ■ , Pk) be the Frobenius notation for a partition. Then, 

Hi^iati - atj)(pi - Pj) 1 

S(a|/3)(too) = — ^k , " , ^ TTT^fc ~) ' ^ 388 ) 



T k r„ „ \la a \ k / x k 



, ^ IK. ,'»• - - ft) A A (-) ft (-q) ft fwn 



Sf Q im(t(oo,q)) = : t r , (390) 

nL=i(g- ft - g Qj+1 ) n- =1 (g;g) Q , n- =1 (g;g) ft 

S( a \(s)(t(a,q)) = (391) 

n^-(g"' +1 - g a ' +1 )(g-* - g- ft ) * A HN^(^ 

n- ; =i(g- ft -g Q3+1 ) ti fi (g;g) ft 



From (|379[) and from (|388[) it may be derived 
Lemma 5. Let A = (Ai, . . . , A&) be a strict partition. Then 



. f < A ) A 1 A Ai — Ai 



i= 1 i<j 
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E.5 "Neutral" two-component BKP hierarchy and 2-BKP hierarchy 

Two-component sBKP. Consider the following tau function of the "small" two-component BKP 
hierarchy 

T (B \t',t",C) = (0,0|r B (t')f B (t")ex P CWn<L|0,0) (394) 

n , m > 

= £ 2" W-i'(«det C X Q a Qp (395) 
Aep ^ ' ^ ' 

where A = (a\f3) and C\ :— det {C) a . p. and Q Q (|) is the projective Schur function related to a (strict) 

partition a, see [TB]- Neutral fermions n ,</>„ were considered in p58[) . If we take C nrn — C nrn (t) = 
e m ~ ™S( n i m )(t) we obtain 

r( B )(t',t",C(t)) = g 2 -tW e u m -u {a} Qa ^ q ^ sx{t) (396) 

Thus this tau function is a tau function of three hierarchies at the same time which are the 1DKP and 
two "small" BKP ones. 

Let us present the following 'symmetric' fermionic representation of this tau function which looks 
little in common with usual notation for tau functions 

(o, o, o|r(t)r s (t')r B (t")e i: ™.'"^ ef "" ff "*»*-.-^»#- |o, o, o) (397) 

2-sBKP tau function. This is 

r(t',i',A) = (0|r B (t')e E »- ^"^f s (t')|0) (398) 
where A is an anti-symmetric matrix. Thanks to 

(0|r B (t')= £ («l 2-*^«> Q Q , r B (F)|0)= \a)i^ m Qp, (a\fl=S~,f> (399) 
oedp /3gdp 

by Wick theorem it may be written in form of double series over strict partitions as 

= £ 2-^-^A a , p Q a (t')Qp(t') (400) 

Q,0GDP 

where 

A a ^ = det [A anJm ] , A nm = (0\(f>- n e^ Ai ^ ( f> m \Q) (401) 

F Appendix. Simplest 1DKP solitons 

We want to present certain types of solitonic solutions typical for 1DKP hierarchy. 

I. KP hierarchy is a reduction of 1DKP one, therefore KP solitons are also 1DKP ones. First let me 
remind a typical KP two-soliton tau function (which is also 1DKP two-soliton tau function which will be 
denoted by r| s0 ,(t)) 

rLol^) = (0|r(t)e aiWpiWt(9l)+Q2 ^ (p2 ^ t(92) |0) = (402) 
1 + e 5l e^-=i -«n + e *" e ££=i + e Sl2 ]J e 5< e^™=i '-W'-O ( 40 3) 

i=l,2 

where 



5i = fog 8 l3 = "r (404) 



Oi s = (Pi -Pj)(Qi ~Qj) 

'Pi-qi' 13 ipi - Qj)(Pj - Qi) 

This two-soliton solution exhibit resonance behavior [43] when p\ — > p2, the same occurs when q\ — > qi 
(such solitons has the Manakov's Y-form). 
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For 1DKP one can present a larger amount of various one- and two-soliton tau functions. Examples 
are written down below. 

II. One of them we may obtain by replacing the right vacuum vector |0) by |f2) in the expectation in 



the right hand side of (|402|) : 

rlLii*) = (0|r(t)e Ql ^ (pi) ^ t( « l)+a2 ^ P2 ^ t(92) |f}) = (405) 

r (t) I 1 + e 5l+Al e^™=i t "»(Pi > -9i*) + e 5 2+ A 2 e E~=i * m (p™-92") _|_ e <5i2+Ai2 J-J e Si+K e T,^tm{pT-iT) 

(406) 

where 5 1,2 and <5i2 are given by (|404[) and where 



This solution describes interaction of two solitons, i — 1,2, each is 

riUt) = r (t) (l + e^e^^^-O) , i = 1,2 (408) 

(notice the To(t) factor in front of the right-hand side). 

Now two-soliton solution exhibit resonance behavior also in case when p\ —> q% , the same occurs 
when P2 — > 1 . 



One can say that we add 'KP solitons' to 1BKP background solution p22[) . 

III. Adding of 'KP solitons' to a 1BKP background solution r (t) of form (|204|) yields 



<w(t) = (N + l\T(t)e a ^^g-\l) = r (t) (l + aj/f'e) (409) 
IV. 1DKP 1-soliton solution may be also as follows 

'"'•"■> u; = 1 + e" lz 

i=l,2 



Tu i(t) = (0\T(t)e ai ^ PlWp2)+a2 ^ (qi) ^ iq2) \0) = 1 + e Sl2 JJ e*e^> ^W*-!!™) (410) 



This soliton is characterizes by four spectral parameters p\,p2,qi 1 <?2 instead of a pair unlike, say, typical 
KP soliton. 

The corresponding '2-soliton' tau function is 



T. 
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(t) = (0|r(t)e ai2 ^' Pl ^^ P2 ^ +o i 2 ^ + ' 9l ^ t ^ 2 ^ +a34 ^^ P3 ^' P4 ^ +a 3 4 ^ t ^ 3 ' 1 '' t ^ 4 - l |0) = (411) 



1 + p/12 , 12 Yl eE^=iPTt m JJ e -E~=i flT*-. +e ^4,34 J| e E~ =1 pr*m J"J e -E~ =1 ?r« m (412) 

i=l,2 i=l,2 i=3,4 i=3,4 

+ g£12 ' 34 1 I e^™^^*"" e-^™^ 8 ™*" + e E34 ' 12 [J e E ~=i^ tm ] [ e - £™=i 9 *' t "» + (413) 

i=l,2 i=3,4 i=3,4 i=l,2 

+ e ei234 TJ e E~ =1 * m (pr-?D (414) 

1=1,2,3,4 

where 

£ij,kl = - log (ft - qk)(pi ~ Ql)(pj - Qk)(Pj - Ql) + Oij + log(ft - pj) + a* kl + \og{q k - qi) (415) 
and where 

nt<j(Pi-Pj)(gi-g;) , 1 , *s ,„ 1R , 

£1234 = =3 ; h log(ai 2 a 12 a 3 4a 34 ) (416) 

rii,,-=i(pj - ft) 

r /y (t) = {Q\T(t)e^™ c ™^ Pn) ^ {qm) p) = ( 417 ) 
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e e n n,t<3 ' fc( r"' Pm,)(g "' qni) ( 4is ) 

k=Omi,...,mi»i,-,»i i<j = l (Pm ( — (Zn.,) 

V. In case N > 0: 
Onc-soliton 

rL;(t) - (7V|r(t) e ^ t fe^ t fe) 3 -|0) (419) 

= ]T s X (t)+a( qi q 2 ) N Yl fi A(t + fa" 1 ] + fe" 1 ]) (420) 

e(\)<N e(\)<N+2 

Now, let us be interested in the case where 



1 W 

t m = Ux) = - E ( 421 ) 



»=i 
Then 
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